Non-Abelian statistics and topological quantum information processing in ID wire networks 
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Topological quantum computation provides an elegant way around decoherence, as one encodes quantum 
information in a non-local fashion that the environment finds difficult to corrupt. Here we establish that one of 
the key operations — braiding of non-Abelian anyons — can be implemented in one-dimensional semiconductor 
wire networks. Previous worlM^ provided a recipe for driving semiconducting wires into a topological phase 
supporting long-sought particles known as Majorana fermions that can store topologically protected quantum 
information. Majorana fermions in this setting can be transported, created, and fused by applying locally tunable 
gates to the wire. More importantly, we show that networks of such wires allow braiding of Majorana fermions 
and that they exhibit non-Abelian statistics like vortices in a p + ip superconductor. We propose experimental 
setups that enable the Majorana fusion rules to be probed, along with networks that allow for efficient exchange 
of arbitrary numbers of Majorana fermions. This work paves a new path forward in topological quantum com- 
putation that benefits from physical transparency and experimental realism. 

PACS numbers: 



The experimental realization of a quantum computer ranks 
among the foremost outstanding problems in condensed mat- 
ter physics, particularly in light of the revolutionary rewards 
the achievement of this goal promises. Typically, decoher- 
ence presents the primary obstacle to fabricating a scalable 
quantum computer. In this regard topological quantum com- 
puting holds considerable promise, as here one embeds quan- 
tum information in a non-local, intrinsically decoherence-free 
fashion^ ^. The core ideas can be illustrated with a toy model 
of a spinless, two-dimensional (2D) p + ip superconductor. 
Vortices in such a state bind exotic particles known as Ma- 
jorana fermions, which cost no energy and therefore gener- 
ate a ground state degeneracy. Because of the Majoranas, 
vortices exhibit non-Abelian braiding statisticJ^lGH adiabat- 
ically exchanging vortices noncommutatively transforms the 
system from one ground state to another. Quantum informa- 
tion encoded in this ground state space can be controllably 
manipulated by braiding operations — something the environ- 
ment finds difficult to achieve. 

Despite this scheme's elegance, realizing suitable topo- 
logical phases poses a serious challenge. Most effort has 
focused on the quantum Hall state at filling fractiorl'^'--' 
ly = 5/2, though very recently the list of candidate ex- 
perimental sy stem s has rapidly expanded. Indeed, topolog- 
ical insulator^'^''^', semiconductor heterostructures'^ non- 
centrosymmetric superconductoriSEl and quantum Hall sys- 
tems at integer plateau transitions^^ can all be engineered 
into non-Abelian topological phases similar to a spinless 
p + ip superconductor. More recently, two groups'^ rec- 
ognized that topological superconductivity can be perhaps 
most easily engineered in one-dimensional (ID) semiconduct- 
ing wires deposited on an s-wave superconductor. These 
proposals provide the first realistic experimental setting for 
Kitaev's^'' ID topological superconducting state. Remarkably, 
the ends of s uch w ires support a localized, zero-energy Majo- 
rana fermiorPS^H Motivated by the exciting possibility of ex- 



perimentally realizing this phase, we examine the prospect of 
exploiting ID semiconducting wires for topological quantum 
computation. 

The suitability of ID wires for this purpose is by no means 
obvious. One of the operations desired for topological quan- 
tum computation is braiding (though measurement-only ap- 
proaches sidestep the need to physically braid particles^), 
and here non-Abelian statistics is critical. While Majorana 
fermions can be transported, created, and fused in a physically 
transparent fashion by applying independently tunable gates 
to the wire, braiding and, in particular, non-Abelian statistics 
poses a serious puzzle. Indeed, conventional wisdom holds 
that braiding statistics is ill-defined in ID, since particles must 
pass through one another at some point during the exchange. 
This problem can be surmounted by considering networks of 
ID wires, the simplest being the T-junction of Fig. [3] Even in 
such networks, however, non-Abelian statistics does not im- 
mediately follow as recognized by Wimmer et al}'^ For ex- 
ample, in a 2D p + ip superconductor, vortices binding the 
Majora nas pl ay an integral role in establishing non-Abelian 
statistic^^El We demonstrate that despite the absence of vor- 
tices in the wires, Majorana fermions in semiconducting wires 
exhibit non-Abelian statistics and transform under exchange 
exactly like vortices in a p + ip superconductor. 

We further propose experimental setups ranging from min- 
imal circuits (involving one wire and a few gates) that can 
probe the Majorana fusion rules, to scalable networks that en- 
able efficient ex change of many Majoranas. The 'fractional 
Josephson effecti^lSEIISl along with Hassler et aUs recent 
proposal^^ enable readout of the topological qubits in this set- 
ting. The relative ease with which ID wires can be driven into 
a topological superconducting state, combined with the phys- 
ical transparency of the manipulations, render the setups dis- 
cussed here extremely promising venues for topological quan- 
tum information processing. Although braiding of Majoranas 
alone does not permit universal quantum computation** ■22j'2ni 
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implementation of the ideas introduced here would constitute 
a critical step towards this ultimate goal. 
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I. MA JORANA FERMIONS IN ' SPINLESS ' p-WAVE 
SUPERCONDUCTING WIRES 

We begin by discussing the physics of a single wire. Valu- 
able intuition can be garnered from Kitaev's toy model for a 
spinless, p-wave superconducting iV-site chain--; 



(b) 



^ B ^ ^ 



semiconducting wire 



N 



N-1 



H = -n^ cic^ - ^ {tclc^+i + I A|e**Ca;C:r+i + h.c.) ^-wave superconductor 



x=l 



x=l 



f 




(1) 

where Cx is a spinless fermion operator and /i, t > 0, and 
|A|e''^ respectively denote the chemical potential, tunneling 
strength, and pairing potential. The bulk- and end-state struc- 
ture becomes particularly transparent in the special casd^ 
fi = 0,t = I A|. Here it is useful to express 
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with ^a.x — lax Majorana fermion operators satisfying 
ila.xila' ,x'} ~ "^Saa'Sxx'- Thcsc cxprcssions cxposc the 
defining characteristics of Majorana fermions — they are their 
own antiparticle and constitute 'half of an ordinary fermion. 
In this limit the Hamiltonian can be written as 
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Consequently, ib.x and ia,x+i combine to form an ordi- 
nary fermion dx = {lA,x+i + i7-B,x)/2 which costs en- 
ergy 2t, reflecting the wire's bulk gap. Conspicuously ab- 
sent from H, however, are ^a,i and 1b,n, which represent 
end-Majorana modes. These can be combined into an ordi- 
nary (though highly non-local) zero-energy fermion dcnd = 
(7a, 1 + *7-B,7v)/2. Thus there are two degenerate ground 
states |0) and |1) — d\^^\Q), where dcndlO) ~ 0, which serve 
as topologically protected qubit states. Figure [TJ a) illustrates 
this physics pictorially. 

Away from this special limit the Majorana end states no 
longer retain this simple form, but survive provided the bulk 
gap remams finitd^Il. This occurs when < 2t, where a 
partially filled band pairs. The bulk gap closes when |/^| = 2t, 
and for larger a topologically trivial superconducting state 
without end Majoranas emerges. Here pairing occurs in either 
a fully occupied or vacant band. 

Realizing Kitaev's topological superconducting state exper- 
imentally requires a system which is effectively spinless — 
i.e., exhibits one set of Fermi points — and p-wave pairs at the 
Fermi energy. Both criteria can be satisfied in a spin-orbit- 
coupled semiconducting wire deposited on an s-wave super- 
conductor by applying a magnetic field' ^ [see Fig.[T];b)]. The 



FIG. 1: (a) Pictorial representation of the ground state of Eq. ([TJ in 
the Umit ^ = Q, t = |A|. Each spinless fermion in the chain is 
decomposed in terms of two Majorana fermions ia,x and yB,x- Ma- 
joranas 7s, I and 7a,i+i combine to form an ordinary, finite energy 
fermion, leaving two zero-energy end Majoranas 7,4,1 and 7b, -iv as 
showiP. (b) A spin-orbit-coupled semiconducting wire deposited on 
an s-wave superconductor can be driven into a topological supercon- 
ducting state exhibitinesuch end Majorana modes by applying an 
external magnetic fiel&' . (c) Band structure of the semiconducting 
wire when B = (dashed lines) and B 7^ (solid lines). When fi 
lies in the band gap generated by the field, pairing inherited from the 
proximate superconductor drives the wire into the topological state. 



simplest Hamiltonian describing such a wire reads 
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The operator 'ijj^x corresponds to electrons with spin a, effec- 
tive mass TO, and chemical potential /i. (We suppress the spin 
indices except in the pairing term.) In the third term, u denotes 
the (Dresselhaus^' and/or Rashba^^) spin-orbit strength, and 
<j = (cr^, CT*', cr^) is a vector of Pauli matrices. This coupling 
favors aligning spins along or against the unit vector e, which 
we assume lies in the {x, y) plane. The fourth term represents 
the Zeeman coupling due to the magnetic field < 0. Note 
that spin-orbit enhancement can lead to^^ g ^ 2. Finally, the 
last term reflects the spin-singlet pairing inherited from the 
s-wave superconductor via the proximity effect. 

To understand the physics of Eq. Q, consider first Bz — 
A = 0. The dashed lines in Fig. ^c) illustrate the band 
structure here — clearly no 'spinless' regime is possible. In- 
troducing a magnetic field generates a band gap cx j-B^j at 
zero momentum as the solid line in Fig. [TJc) depicts. When 
/i lies inside of this gap the system exhibits only a single pair 
of Fermi points as desired. Turning on A which is weak com- 
pared to the gap then effectively p-wave pairs fermions in the 
lower band with momentum k and — fc, driving the wire into 
Kitaev's topological phase' ^. [Singlet pairing in Eq. ^ gen- 
erates p-wave pairing because spin-orbit coupling favors op- 
posite spins for k and — fc states in the lower band.] Quan- 
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titatively, realizing the topological phase require^**^ |A| < 
g^B\Bz\/'2,, which we hereafter assume holds. The opposite 
limit |A| > effectively violates the 'spinless' cri- 

terion since pairing strongly intermixes states from the upper 
band, producing an ordinary superconductor without Majo- 
rana modes. 

In the topological phase, the connection to Eq. ([T]l becomes 
more explicit when g^B \Bz\ ^ "rnu^ , | A | where the spins 
nearly polarize. One can then project Eq. Q onto a sim- 
pler one-band problem by writing ijj^x ^ "^^^^^J^^^^^a; '^^^ 
i^ix ^ with '5 J, the lower-band fermion operator. To 
leading order, one obtains 
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where /ios 
field reads 



M + .9/^51-8^1/2 and the effective p-wave pair 
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The dependence of ip^s on e will be important below when 
we consider networks of wires. Equation (j5]l constitutes an ef- 
fective low-energy Hamiltonian for Kitaev's model in Eq. ([T]l 
in the low-density limit. From this perspective, the existence 
of end-Majoranas in the semiconducting wire becomes mani- 
fest. We exploit this correspondence below when addressing 
universal properties such as braiding statistics, which must be 
shared by the topological phases described by Eq. Q and the 
simpler lattice model, Eq. ([T]l. 

We now seek a practical method to manipulate Majorana 
fermions in the wire. As motivation, consider applying a gate 
voltage to adjust /i uniformly across the wire. The excitation 
gap obtained from Eq. Q at fc = varies with fi via 



£^gap(fc = 0) = 



g^J■B\Bz 
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For l/^l < /ic = yJ{g^BBz/'i)'^ — |Ap the topological phase 
with end Majoranas emerges, while for > a topologi- 
cally trivial phase appears. Applying a gate voltage uniformly 
thus allows one to create or remove the Majorana fermions. 
However, when = /ic the bulk gap closes, and the exci- 
tation spectrum at small momentum behaves as i?gap(fc) ~ 
/it;|fc|, with velocity v = 2u\/S.\/ {giiB\Bz\)- The gap closure 
is clearly undesirable, since we would like to manipulate Ma- 
jorana fermions without generating additional quasiparticles. 

This problem can be circumvented by employing a 'key- 
board' of locally tunable gates as shown in Fig. |2] each 
of which impacts over a finite length igatc of the wire. 
When a given gate locally tunes the chemical potential across 
— lie, a finite excitation gap i?gap ^ ^wTr/igatc remains. 
(Roughly, the gate creates a potential well that supports only 
k larger than ^ 7r/Lgato-) Assuming giJ,B\Bz\/2 ^ 2|A| and 
hu ~ O.leVA yields a velocity v ^ lO'^m/s; the gap for a 
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FIG. 2: Applying a 'keyboard' of individually tunable gates to the 
wire allows one to locally control which regions are topological (dark 
blue) and non-topological (light blue), and hence manipulate Majo- 
rana fermions while maintaining the bulk gap. As (a) and (b) illus- 
trate, sequentially applying the leftmost gates drives the left end of 
the wire non-topological, thereby transporting 71 rightward. Nucle- 
ating a topological section of the wire from an ordinary region or vice 
versa creates pairs of Majorana fermions out of the vacuum as in (c). 
Similarly, removing a topological region entirely or connecting two 
topological regions as sketched in (d) fuses Majorana fermions into 
either the vacuum or a finite-energy quasiparticle. 



0.1 /im wide gate is then of order IK. We consider this a con- 
servative estimate — heavy -element wires such as InSb and/or 
narrower gates could generate substantially larger gaps. 

Local gates allow Majorana fermions to be transported, cre- 
ated, and fused as outlined in Fig. [2] As one germinates pairs 
of Majorana fermions, the ground state degeneracy increases 
as does our capacity to topologically store quantum informa- 
tion in the wire. Specifically, 2n Majoranas generate n ordi- 
nary zero-energy fermions whose occupation numbers specify 
topological qubit states. Adiabatically braiding the Majorana 
fermions would enable manipulation of the qubits, but is not 
possible in a single wire. Thus we now turn to the simplest ar- 
rangement which allows for exchange — the T-junction of Fig. 

El 



II. MAJORANA BRAIDING AND NON-ABELIAN 
STATISTICS 

First, we explore the physics at the junction where the wires 
in Fig. |3] meet (see the Supplementary Material for a more 
detailed discussion). It will be useful to view the T-junction 
as composed of three segments whose ends meet at a point. 
When only one segment realizes a topological phase, a single 
zero-energy Majorana fermion exists at the junction. When 
two topological segments meet at the junction, as in Figs.[3|a) 
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FIG. 3: A T-junction allows for adiabatic exchange of two Majorana 
fermions bridged by either a topological region (dark blue lines) as 
in (a)-(d), or a non-topological region (light blue lines) as in (e)-(h). 
Transport of Majorana fermions is achieved by gates as outlined in 
Fig. [2] The arrows along the topological regions in (a)-(d) are useful 
for understanding the non-Abelian statistics as outlined in the main 
text. 



and (b), generically no Majorana modes exist there. To see 
this, imagine decoupling the two topological segments so that 
two Majorana modes in close proximity exist at the junction; 
restoring the coupling generically combines these Majoranas 
into an ordinary, finite-energy fermion. 

As an illustrative example, consider the setup of Fig. [3]^a) 
and model the left and right topological segments by Kitaev's 
model with ^ = and t — |A| in Eq. ([T]). [For simplic- 
ity we will exclude the non-topological vertical wire in Fig. 

ta).] Suppose furthermore that the cj) — (p / ^ in the left/right 
ains and that the fermion cl,n at site of the left chain 
couples weakly to the fermion cr^i at site 1 of the right chain 
via Hy = — r(c^ n'^R,! + h.c). Using Eq. (jij, the end Ma- 



joranas at the junction couple as follows 
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and therefore generally combine into an ordinary fermion^^. 
An exception occurs when the regions form a 7r-junction — 
that is, when (t>L — 4'r = tt — which fine-tunes their coupling 
to zero. Importantly, coupling between end Majoranas in the 



semiconductor context is governed by the same (f>L — 4>r de- 
pendence as in Eq. (jsji' ^. 

Finally, when all three segments are topological, again only 
a single Majorana mode exists at the junction without fine- 
tuning. Three Majorana modes appear only when all pairs of 
wires simultaneously form mutual tt junctions (which is pos- 
sible as described in the Supplementary Material, since the 
superconducting phases are defined with respect to a direction 
in each wire). Recall from Eq. (|6]l that the spin-orientation 
favored by spin-orbit coupling determines the effective super- 
conducting phase of the semiconducting wires. Two wires 
at right angles to one another therefore exhibit a 7r/2 phase 
difference, well away from the pathological limits mentioned 
above. 

The T-junction permits exchange of Majoranas from either 
the same or different topological regions. First, consider the 
configuration of Fig.|3ja) where the horizontal wire resides in 
a topological phase while the vertical wire is non-topological. 
Counterclockwise exchange of 71 and 72 can be implemented 
as outlined in Figs. [3jb)-(d). Here, one shuttles 71 to the 
junction by making the left end non-topological; transports 
7i downward by driving the vertical wire into a topological 
phase; transports 72 leftward in a similar fashion; and finally 
directs 71 up and to the right. Exchange of two Majorana 
fermions connected by a non-topological region as in Fig.[3|e) 
can be similarly achieved — counterclockwise exchange of 71 
and 72 proceeds as sketched in Figs.|3jf)-(h). 

While the Majoranas can now be exchanged, non-Abelian 
statistics is not obvious in this context. Recall how non- 
Abelian statist ics o f vortices arises in a spinless 2D p + ip 
superconductojE^IiH following Ivanov's approach. Ultimately, 
this can be deduced by considering two vortices which bind 
Majorana fermions 71 and 72. Since the spinless fermion 
operators effectively change sign upon advancing the super- 
conducting phase by 27r, one introduces branch cuts emanat- 
ing from the vortices; crucially, a Majorana fermion changes 
sign whenever crossing such a cut. Upon exchanging the 
vortices, 72 (say) crosses the branch cut emanating from the 
other vortex, leading to the transformation rule 71 72 
and 72 —71 which is generated by the unitary operator 
U12 — exp(7r727i/4). With many vortices, the analogous 
unitary operators Uij corresponding to the exchange of 7^ and 
7j do not generally commute, implying non-Abelian statistics. 

Following an approach similar to that of Stern et al?'^, we 
now argue that Majorana fermions in semiconducting wires 
transform exactly like those bound to vortices under exchange, 
and hence also exhibit non-Abelian statistics. This can be 
established most simply by considering the exchange of two 
Majorana fermions 71 and 72 as illustrated in Figs. |3ja)-(d). 
At each step of the exchange, there are two degenerate ground 
states |0) and |1) = /tlO), where / = (71 + ^72)72 annihi- 
lates |0). In principle, one can deduce the transformation rule 
from the Berry phases Xn = Im / dt{n\dt\n) acquired by the 
many-body ground states \n) = |0) and though in practice 
these are hard to evaluate. 

Since exchange statistics is a universal property, however, 
we are free to deform the problem to our convenience pro- 
vided the energy gap remains finite. As a first simplification. 



since the semiconductor Hamihonian and Kitaev's model in 
Eq. ([T]l can be smoothly connected, let us consider the case 
where each wire in the T-junction is described by the latter. 
More importantly, we further deform Kitaev's Hamiltonian to 
be purely real as we exchange 71.2- The states |0) and |1) 
can then also be chosen real, leading to an enormous simpli- 
fication: while these states still evolve nontrivially the Berry 
phase accumulated during this evolution vanishes. 

For concreteness, we deform the Hamiltonian such that 
/i < and i = A = in the non-topological regions of 
Fig. |3] For the topological segments, reality implies that the 
superconducting phases must be either or tt. It is useful to 
visualize the sign choice for the pairing with arrows as in Fig. 
[3] (To be concrete, we take the pairing | A|e*'^CjCj+i such 
that the site indices increase moving rightward/upward in the 
horizontal/vertical wires; the case = then corresponds 
to rightward/upward arrows, while leftward/downward arrows 
indicate ~ tt.) To avoid generating tt junctions, when two 
topological segments meet at the junction, one arrow must 
point into the junction while the other must point out. With 
this simple rule in mind, we see in Fig.|3]that although we can 
successfully swap the Majoranas while keeping the Hamilto- 
nian real, we inevitably end up reversing the arrows along the 
topological region. In other words, the sign of the pairing has 
flipped relative to our initial Hamiltonian. 

To complete the exchange then, we must perform a gauge 
transformation which restores the Hamiltonian to its original 
form. This can be accomplished by multiplying all fermion 
creation operators by i; in particular, = (71 — 1^2) I '2, -> 

~ (72 + *7i)/2. It follows that 71 — > 72 and 72 
which the unitary transformation U12 — exp(7r727i/4) gen- 
erates as in the 2D p + ip case. We stress that this result 
applies also in the physically relevant case where gates trans- 
port the Majoranas while the superconducting phases remain 
fixed; we have merely used our freedom to deform the Hamil- 
tonian to expose the answer with minimal formalism. Ad- 
ditionally, since Figs. [3|e)-(h) also represent a counterclock- 
wise exchange, it is natural to expect the same result for this 
case. The Supplementary Material analyzes both types of ex- 
changes from a complementary perspective (and when the su- 
perconducting phases are held fixed), confirming their equiva- 
lence. There we also establish rigorously that in networks sup- 
porting arbitrarily many Majoranas exchange is implemented 
by a set of unitary operators Uij analogous to those in a 2D 
p + ip superconductor Thus the statistics is non-Abelian as 
advertised. 



III. DISCUSSION 

The keyboard of gates shown in Fig. [2] and the T-junction 
of Fig. [3] provide the basic elements allowing manipulation 
of topological qubits in semiconducting wires. In principle, a 
single T-junction can support numerous well-separated Majo- 
rana modes, each of which can be exchanged with any other 
(First, create many Majoranas in the horizontal wire of the 
T-junction. To exchange a given pair, shuttle all intervening 
Majoranas down to the end of the vertical wire and then carry 
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FIG. 4: (a) Example of a semiconductor wire network which al- 
lows for efficient exchange of many Majorana fermions. Adjacent 
Majoranas can be exchanged as in Fig. [3] while non-adjacent Majo- 
ranas can be transported to the lower wire to be similarly exchanged, 
(b) Minimal setup designed to detect the non-trivial Majorana fusion 
rules. Majoranas 71,2 are first created out of the vacuum. In the left 
path, 72 is shuttled rightward, and Majoranas 73,4 always combine 
to form a finite-energy state which is unoccupied. In the right path, 
73,4 are also created out of the vacuum, and then 72 and 73 fuse with 
50% probability into either the vacuum or a finite-energy quasiparti- 
cle. The Josephson current flowing across the junction allows one to 
deduce the presence or absence of this extra quasiparticle. 



out the exchange using the methods of Fig. [3]) However, net- 
works consisting of several T-junctions — such as the setup of 
Fig. |4|a) — enable more efficient Majorana exchange. In the 
figure, all adjacent Majorana fermions can be immediately 
swapped using Fig.|3] while non-adjacent Majoranas can be 
shuttled down to the lower wire to be exchanged. This 'lad- 
der' configuration straightforwardly scales up by introducing 
additional 'rungs' and/or 'legs'. 

As Fu and Kane suggested in the topological insulator 
contexfiS fusing Majorana fermions across a Josephson junc- 
tion provides a readout method for the topological qubit states. 
We illustrate the physics with the schematic setup of Fig.|4|b), 
which extends the experiments proposed in Refs. lII2l to allow 
the Majorana fusion rules to be directly probed. Here a semi- 
conducting wire bridges two s-wave superconductors with ini- 
tial phases tpl/ij; we assume A</5* = ip\ — ip]-^ 7^ tt. Three 
gates drive the wire from an initially non-topological ground 
state into a topological phase. Importantly, the order in which 
one applies these gates qualitatively affects the physics. As 
we now discuss, only in the left path of Fig.|4]^b) can the qubit 



6 



state at the junction be determined in a single measurement. 

Consider first germinating Majorana fermions 71 and 72 by 
applying the left gate. Assuming /a = (71 + 172 )/2 initially 
costs finite energy as 71 and 72 separate, the system initializes 
into a ground state with /a unoccupied. Applying the central 
and then right gates shuttles 72 to the other end [see the left 
path of Fig.|4|b)]. Since a narrow insulating barrier separates 
the superconductors, an ordinary fermion Jb = (73 + j74)/2 
arises from two coupled Majoranas 73 4 at the junction. Simi- 
lar to Eq. (jsj, the energy of this mode is well-captured bjPEEIl 
Hj ~ ze'7374 = ^'i^flfB - 1) where = 5cos(Ay)V2) 
with non-universal S. The system has been prepared in a 
ground state, so the fs fermion will be absent if e* > but 
occupied otherwise. 

Suppose we now vary the phase difference across the junc- 
tion away from its initial value to A(p. The measured Joseph- 
son current (see Supplementary Material for a pedagogical 
derivation) will then be 

/ = I ^ = jsgn{e^) sin(A^/2) + h,, (9) 

where E is the ground-state energy and l2e denotes the usual 
Cooper-pair-tunneling contribution. The first term on the right 
reflects single-electron tunneling originating from the Majo- 
ranas 73 4. This 'fractional' Josephson current exhibits An 
periodicity in Aip, but 2tt periodicity in the initial phase dif- 
ference A(p'. 

The right path in Fig.[4]^b) yields very different results, re- 
flecting the nontrivial Majorana fusion rules. Here, after cre- 
ating 7i 2 one applies the rightmost gate to nucleate another 
pair 73.4. Assuming /a and Jb defined as above initially 
cost finite energy, the system initializes into the ground state 
|0,0) satisfying /a/b|0,0) = 0. Applying the centi-al gate 
then fuses 72 and 73 at the junction. To understand the out- 
come, it is useful to re-express the ground state in terms of 
Ia = (71 + i74)/2 and fg = (72 + ns)/'^- In this ba- 
sis |0,0) = (|0',0') - i\l',l'))/V2, where f'^ g annihilate 

|0',0') and |1',1') = /^Vb|0',0'). Following our previous 
discussion, f'g acquires finite energy at the junction, lifting 
the degeneracy between |0',0') and Measuring the 

Josephson current then collapses the wavefunction with 50% 
probability onto either the ground state, or an excited state 
with an extra quasiparticle localized at the junction. In the 
former case Eq. (|9]l again describes the current, while in the 
latter the fractional contribution simply changes sign. 

In more complex networks such as that of Fig.|4|a), fusing 
the Majoranas across a Josephson junction — and in particu- 
lar measuring the sign of the fractional Josephson current — 
similarly allows qubit readout. Alternatively, the interesting 
recent proposal of Hassler et al'^ior reading qubit states via 
ancillary non-topological flux qubits can be adapted to these 
setups (and indeed was originally discussed in terms of an iso- 
lated semiconducting wire^^). 

To conclude, we have introduced a surprising new venue 
for braiding, non-Abelian statistics, and topological quantum 
information processing — networks of one-dimensional semi- 
conducting wires. From a fundamental standpoint, the abil- 



ity to realize non-Abelian statistics in this setting is remark- 
able. Perhaps even more appealing, however, are our pro- 
posal's physical transparency and experimental promise, par- 
ticularly given the feats already achieved in Ref. 35 While 
topological quantum information processing in wire networks 
requires much experimental progress, observing the distinct 
fusion channels characteristic of the two paths of Fig. |4|b) 
would provide a remarkable step en route to this goal. And 
ultimately, if braiding in this setting can be supplemented by 
a 7r/8 phase gate and topological charge measurement of four 
Majoranas, wire networks may provide a feasible path to uni- 
versal quantum computatio na l ^^" 1 



IV. SUPPLEMENTARY MATERIAL 
A. Properties of the T-junction 

Here we investigate in greater detail the properties of the 
junction in Fig. [3] where the three wire segments meet. There 
are three cases to consider, corresponding to the situations 
where one, two, or all three of the wire segments emanat- 
ing from the junction reside in a topological superconduct- 
ing state. It is conceptually simplest to address each case by 
viewing the T-junction as composed of three independent wire 
segments as in Fig. [5] which initially decouple from one an- 
other In this limit a single Majorana exists at the end of each 
topological segment. One can then straightforwardly couple 
the wire segments at the junction and explore the fate of the 
Majorana end states. 

Suppose that the phases of the p-wave pair fields in each 
region are (j)A/B /c shown in Fig. |5] [in the semiconduc- 
tor wire context, these phases correspond to cpcs in Eq. (6) of 
the main text]. To be precise, if the wires are described by 
a lattice model, we define these phases relative to a pairing 
term | A|e*'^°CjCj+i such that the site indices increase moving 
rightward in the horizontal wires and upward in the vertical 
wire. A similar convention can be employed in the semicon- 
ductor wire context. Now suppose we allow single-electron 
tunneling between the ends of each segment, with amplitude 
r as shown schematically in Fig. [5] (Pairing between elec- 
trons residing at the ends of each region is also generally al- 
lowed, but does not change any qualitative results below and 
will therefore be neglected.) For convenience we will assume 
that the tunneling strength is weak compared to the bulk gaps 
in the wires, which will allow us to focus solely on the Ma- 
jorana end states; our conclusions, however, are more general 
and do not require this assumption. 

In the setup of Fig. [Sj^a) with only one topological region, 
the Majorana 71 is qualitatively unaffected by the coupling 
to the non-topological wires. At most its wavefunction can 
be quantitatively modified, but it necessarily remains at zero 
energy. This reflects the familiar topological protection of a 
single isolated Majorana mode in a gapped system. 

With two of the three wires topological as in Fig. |5jb), 
the end Majoranas 71 and 72 generally combine into an or- 
dinary finite-energy fermion, except with fine-tuning. To 
a good approximation, the Majoranas couple through a 
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7l r 



(a) 



Equation ( 10 1 demonstrates that 71 and 72 remain zero-energy 
modes only when the topological wires form a tt junction, i.e., 
when <j>A = 4*3 + TJ"- 

Finally, consider the case shown in Fig.|5]^c) where all three 
segments are topological. Here the Majoranas 71^2,3 couple 



H123 oc -iV 




7l r 72 



(c) 




FIG. 5: T-junction viewed as three wire segments with p- wave super- 
conducting phases <I>a,b.c- The ends of each segment are coupled 
via tunneling with amplitude F as shown, (a) When only one seg- 
ment is topological, the tunneling can not destroy the Majorana 71 
at the junction, (b) Two topological regions meeting at the junction 
leads to the end Majoranas 71 and 72 generally combining into an 
ordinary, finite-energy fermion, unless the topological wires form a 
TT junction, (c) When all the three wires are topological, the Majo- 
ranas 71,2,3 generally combine to form a finite-energy fermion and 
a single topologically protected Majorana. All three Majoranas re- 
main at zero energy only when all three wire segments form mutual 
TT junctions. 



HamiltoniarpEl 



H12 oc —iT cos 



7i72- 



(10) 



This was discussed in the main text in the context of two wires 
described by Kitaev's toy model in a particular limit, but is 
qualitatively rather general — the in periodicity in (f) a — (f>B 
has a topological origin^^. For instance, end Majoranas in 
two topological semiconducting wires coupled through an or- 
dinary region exhibit the same phase dependence as above^. 



, .'A-<PB , 
COS I ) 7i72 



COS I I 7372 + sm I ) 7173 



(11) 



Note the sine function determirung the coupling between 71 
and 73, which arises because of the conventions we chose for 
defining (p^ above. Recall that to make the problem well- 
defined, we needed to define the phases with respect to a par- 
ticular direction in each wire; otherwise there is an ambigu- 
ity of vr in the definition, since for instance \ A\e^'^CjCj+i — 
I Aje^'-'^'^'^-'cj+iCj. We defined the phases such that the site in- 
dices increase upon moving rightward or upward in the wires. 
But this implies that the site indices in both the left and bot- 
tom wires increase upon moving towards the junction, in con- 
trast to all other pairs of wires. It follows that the splitting 
of 71 and 73 is proportional to cos[{(f)A ^ <t>c — ''^)/'A = 
sm[{(f)A — (t>c)l'A- Hence with our conventions a tt junction 
between these two regions actually corresponds to the case 

(t'A = 4>C- 

The Hamiltonian i/i23 implies that 71.2,3 all remain zero- 
energy modes only when (pA — 4'c = + t^, where 
all pairs of wires form mutual tt junctions. (This remains 
true even when coupling to the ordinary gapped states is in- 
cluded.) Aside from this fine-tuned limit, however, H123 al- 
ways supports one zero-energy Majorana mode and one or- 
dinary finite-energy fermion. As an illustration, consider the 
case 0^ = + TT and (pc ^ (pA, so that only the horizontal 
wires form a tt junction. Here the Hamiltonian simplifies to 



H123 oc — ircos 



(j>c -~ 4>i 



73(72 - 7i)- 



(12) 



It follows that the linear combination (71 + 72) /a/2 remains 
a zero-energy Majorana mode, while 73 and (71 — 72)/\/2 
combine into a finite-energy fermion. While here the zero- 
energy Majorana carries weight only on the horizontal wires 
which formed the tt junction, in general its wavefunction will 
have weight on all three segments. 

As we braid Majorana fermions using the methods de- 
scribed in the main text, it is imperative that we avoid gen- 
erating spurious zero-modes at the T-junction. The above dis- 
cussion implies that we are safe in this regard so long as we 
avoid TT junctions. Fortunately, the semiconducting wires we 
considered naturally avoid such situations, since two wires at 
right angles to one another exhibit effective p-wave phases 
that differ by 7r/2 as discussed in the main text. 
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B. Wavefunction approach to Majorana fermion exchange 

In this section we explore in much greater detail the ex- 
change of Majorana fermions in a ID wire networks. We 
once again emphasize the nontrivial nature of the problem 
since non-Abelian statistics in a 2D p + ip superconductor 
is typically understood as arising because of superconduct- 
ing vortices. One might worry that the Majoranas in the 
wires perhaps bind vortices in the neighboring parent s-wave 
superconductor, but this is certainly not the case. This be- 
comes apparent when one recalls how the effective super- 
conducting Hamiltonian for the wire is derived (see, e.g., 
Ref. i36i ). Namely, one considers a Hamiltonian of the form 
H = ifwiro + Hsc + Ht, where ffwirc and iJsc describe the 
wire and superconductor in isolation, and Ht encodes single- 
electron tunneling between the two. Upon integrating out the 
gapped superconducting degrees of freedom assuming a uni- 
form pair field A, one arrives at an effective Hamiltonian for 
the wire which includes proximity induced pairing terms. Any 
phase variations in the parent superconductor's order param- 
eter are ruled out by assumption, yet Majoranas can never- 
theless exist in the wires. Thus developing a physical picture 
for the exchange in this setting poses an extremely important 
issue. 

Our aim here is to provide greater rigor and a complemen- 
tary picture for the discussion presented in the main text for 
how the Majoranas transform under exchange. We will begin 
by constructing the many-body wavefunctions for a general 
ID wire network supporting an arbitrary number of Majorana 
fermions. We will then establish some important general re- 
sults for braiding that rely on minimal assumptions about the 
underlying Hamiltonian. Here there will be some overlap with 
the approach followed by Bonderson et al^, who recently 
revisited the issue of non-Abelian statistics in the fractional 
quantum Hall context. Following this general analysis, we 
will consider again the exchange of Figs.[3];a)-(d) and explore 
how the Majoranas transform when the braid is implemented 
by keeping the superconducting phases fixed, as would be the 
case in practice. Next, we will turn to an analysis of the ex- 
change outlined in Figs.[3je)-(h) and show that, as claimed in 
the main text, this braid transforms the Majoranas in an identi- 
cal fashion to the braid of Figs.[3]^a)-(d). Finally, we will con- 
sider some special examples where the full many-body wave- 
functions can be analyzed during the exchange, allowing us to 
explore important issues such as the overall phase acquired by 
the ground states under braiding. 



1. Construction of degenerate ground state wavefunctions 

Consider a ID wire network with 2 A/ well-separated, lo- 
calized Majorana modes corresponding to operators 7i,....2j\/ 
that satisfy 7^ = 7^ and 7| = 1. By 'well-separated', we 
mean that different Majorana wavefunctions overlap negligi- 
bly with each other Suppose moreover that the pairs 72^-1 
and 72j were germinated from the vacuum. (For example, one 
could start with a non-topological network, generate 71 and 72 
by nucleating a single topological region, then create 73 and 



74 by forming another far-away topological region, etc.) One 
can construct M fermion operators from these via 



(13) 



which correspond to zero-energy modes (up to corrections 
which are exponentially small in the separation between Ma- 
joranas). These modes give rise to 2*^ degenerate ground 
states which can be labeled by the occupation numbers rij — 
0, 1 for the fj fermions. We would like to construct these de- 
generate ground state wavefunctions and understand the ex- 
change of Majorana fermions in ID wire networks from this 
perspective. 

Let us denote the positive-energy Bogoliubov-de Gennes 
quasiparticle operators by da, each of which must annihilate 
the ground states. As usual, the explicit form of these ground 
states is nontrivial because both the Majorana operators and 
da represent linear combinations of the original fermion cre- 
ation and annihilation operators for the ID wire network. By 
construction the wavefunction 



1 



(14) 



with I vac) the vacuum of the original fermion operators and 
J\f the normalization, must constitute one of the degenerate 
ground states since any da clearly annihilates this state. Be- 
cause we pulled 72j-i and 72^ out of the vacuum, we are guar- 
anteed that will be an eigenvector of /J/, with eigenvalue 

How we label this state in terms of the rij occupation num- 
bers is a matter of convention because there is a sign ambigu- 
ity in the Majorana wavefunctions. Specifically, if (pj is the 
Majorana wavefunction corresponding to jj, then — 0^ also 
denotes a legitimate Majorana wavefunction that preserves the 
relations 7j — and 7| = 1. As an example, suppose that 
upon making some particular overall sign choices for 0i.2 we 
find that /i j?/;) — 0; we would then identify with rii — 0. 
Had we chosen the opposite sign for 02, however, we would 
find instead that fllip) =0 and identify this state with ni = 1 
(sending 02 — )■ —02 is equivalent to sending 72 — )■ —72 and 
/i — > fl). We will assume for concreteness that the signs of 
the Majorana wavefunctions have been chosen such that 
corresponds to the ground state with nj — 1, i.e.. 



(15) 



The ground state with nj — can then be written 

I00---0) =/i/2---/m|V'), (16) 

which is manifestly annihilated by all fj. All other ground 
states can be obtained by applying creation operators fj to 
|00 • • • 0) or, equivalently, annihilation operators fj to the state 

Before moving on, we note that the way we combined Ma- 
joranas in Eq. ( 13 1 to construct the ordinary fj fermion oper- 



ators was convenient because of how we assumed the Majo- 
ranas were germinated, but is by no means unique. One could 
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always choose to combine pairs of Majoranas differently and 
construct operators /j and ground states |n']^n2 ■ ■ ■ n^/ ). The 
ground states in this representation would then be related to 
the ground states we defined above simply by a change of ba- 
sis. 



2. General results for Majorana exchange 

We now proceed to obtain some important generic results 
for Majorana exchange that rely on very minimal assump- 
tions about the underlying Hamiltonian for the ID wire net- 
work. Suppose that we wish to exchange 71 and 72. [Below 
it will prove extremely convenient to work in a basis where 
the Majoranas we braid combine into an ordinary fermion. 
This is the case for the basis we introduced above, since 
/i = (71 + *72)/2. If instead we wanted to exchange, say, 
7i and 73, one would first want to change basis and write the 
ground states in terms of operators such as f[ = (71 + i^i-i)/^ 
and /2 = (72 + 174)72, then proceed as we outline below.] 
Let A be the parameter in the Hamiltonian that varies to im- 
plement the exchange of 71 and 72. If A varies from Ainit 
to Afinai during the course of the exchange, then we require 
that i/(Afinai) = ^^(Ainit) SO that the Hamiltonian returns 
to its original form after the braid. (If the Hamiltonian does 
not return to its original form, then we can not make rigorous 
statements about how the wavefunctions transform under the 
exchange.) 

There are two important contributions that one must under- 
stand to analyze the exchange. The first are the Berry phases 
acquked by the degenerate ground states, which follow from 

(nin2 • • • nM\d\\mim2 ■ ■ ■ ium)- (17) 

We have included the possibility that non-trivial off-diagonal 
Berry phases occur, since a priori these need not vanish. If 
we assume that the degenerate ground states return to their 
precise original form at A/, the Berry phases encode all infor- 
mation about the exchange. However, if we relax this assump- 
tion (which will indeed be useful below), then we additionally 
need to compare the explicit changes between the initial and 
final states. In general, the final ground states could represent 
a nontrivial linear combination of the initial ground states, so 
we write 

\nin2---nM)f = ^ {nj})|mim2 • • • mM)i, 

{rrij} 

(18) 

where the subscripts i/f denote the initial/final states. 
The combination of the Berry phases and the coefficients 
Pll'^^j}; {%}) fully specify the outcome of the exchange. 

While the problem appears daunting, a remarkable amount 
of progress can in fact be made on very general grounds. We 
will make only one additional assumption about the physi- 
cal system. Specifically, to carry out the exchange we will 
assume that only local terms in the Hamiltonian — such as lo- 
cal chemical potentials — need to be modified, and that such 



modifications only impact the Majoranas 71.2 which are be- 
ing braided. This will certainly be the case, for example, in the 
eight-Majorana configuration from Fig. 4(a) of the main text; 
in fact, there any pair can be exchanged without disturbing 
any other Majoranas in the system. However, one is not guar- 
anteed that this is always immediately possible, since there 
may be intervening Majoranas that prevent 71 and 72 from 
being so exchanged. Such cases can be treated in one of two 
ways. First, one can always first transport the Majoranas in the 
system (but importantly without braiding any of them) to pro- 
duce an arrangement where 71 2 can be directly exchanged. 
Alternatively, one can always pairwise braid Majoranas in the 
system until the intervening Majoranas are removed. We will 
proceed below by assuming that in such cases, the prerequisite 
translations or exchanges have already been carried out. 

We will first demonstrate that all off -diagonal Berry phase 
elements vanish trivially in the basis we have chosen, and that 
the exchange depends only on the occupation number ni for 
the /i = (71 +«72)/2 fermion. This follows from our assump- 
tion above that the exchange can be implemented by simply 
adjusting local terms in the Hamiltonian that affect only the 
regions near 71 2. Assuming all other Majoranas are well- 
localized far from 71.2, their wavefunctions will be unaffected 
by these local changes to the Hamiltonian (up to exponentially 
small corrections). Consequently, to an excellent approxima- 
tion f2....,M can be assumed A-independent. 

In contrast, /i and the bulk quasiparticle operators da do 
depend on A. Explicitly displaying the A-dependence, the 
ground state |00 • • • 0) can then be written as 

I00---0) =/i(A)/2---/m|V'(A)). (19) 

Again, all other ground states can be obtained from this by 
application of /i(A)^ and Because /2,....m are A- 

independent, differentiating the states \n1n2 ■ ■ ■ um) with re- 
spect to A does not change the occupation numbers for the 
corresponding zero-energy modes. It therefore follows that 

(711^2 • • • nM\d\\mim2 ■ ■ ■ rriAi) oc (5„2,,„2 • • • (5„J^,^„„ .(20) 

One can see this relation formally by inserting the ground 
states obtained from Eq. ( [T9] l into the left-hand side of the 
above equation. We furthermore have 

{nin2 ■ ■ ■nM\d\\min2 ■ ■■um) oc S^^^rnt (21) 

since with rii ^ mi the above bra and ket are trivially orthog- 
onal because they exhibit different fermion parity. Thus all 
off-diagonal Berry phases indeed vanish, and the only possi- 
bly non-zero Berry phases follow from 

Ao(A) = {0n2- ■ ■nM\d\\Qn2- ■ -um) 

= (VXA)|/i(A)t9;,/i(A)|V'(A)), (22) 

Ai(A) = {lii2- ■ ■nM\dx\ln2- ■ -um) 

= {i^{\)\dxm\)), (23) 

which manifestly depend only on ni. As one would intuitively 
expect, the additional far-away Majorana modes do not impact 
the outcome of the exchange of 712. These modes trivially 
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factor out in the above sense, and the exchange simply sends 
ji -> 7i for i 7^ 1 or 2. 

We can actually evaluate the Berry phases even further. 
From Eqs. ([22]l and ([23]l, we obtain the following relation: 



Ao(A) = Ai(A) + (^(A)|/i(A)naA/i(A)]|VXA)), (24) 

where on the right side the derivative acts only on /i (A). One 
can always decompose d\fi{X) as 

d^hiX) = A(A)/i(A) + A'(A)/i(A)t 

+ ^ [S„(A)d„(A) + s;,(A)d„(A)t] . (25) 

a 

(The right-hand side does not involve f2....,M because by as- 
sumption the corresponding wavefunctions have negligible 
weight near 71 and 72.) Inserting this expression into Eq. ( [24| , 
one finds that only the A{X) term contributes: 



Ao(A) = Ai(A) + A(A). 



(26) 



It is useful to isolate A{X) in Eq. (25 1 by considering the fol- 
lowing anticommutator. 



{/i(A)^9a./i(A)} = A(A). 



(27) 



Suppose now that throughout the exchange the following stan- 
dard relations hold 



/i(A) 



1 



[7i(A) + «72(A)] 



7, (A) - 7j(A)^ 
7,(A)2 - 1. 



(28) 



Note that this does not constitute an additional physical as- 
sumption, but rather a convention; one can always choose the 



operators in this way. The anticommutator in Eq. ( 27 1 is then 



{/i(A)^9a/i(A)} = \[{li,dxji} + {^2,8x12} 

+ i{li,dxj2} ~ i{l2,dxji}],(29) 

where for brevity on the right-hand side we suppressed the A 
dependence. The first two terms vanish trivially because 



{lj,d\lj} = Ijidx-fj) + {dx-fjh] 
= dxh.j,) ^ dxil) ^ 0. 



(30) 



Interestingly, the last two terms in Eq. ( 29 1 also vanish (up to 



exponentially small corrections) because 71 and 72 are local- 
ized far apart from one another. (To see this explicitly, one 
can expand 71 ^ in terms of the original fermion creation and 
annihilation operators for the wire network, evaluate 9a 71, 2 
in terms of this expansion, and then evaluate the anticommu- 
tators.) Thus we obtain the remarkable result that as long as 



Eqs. (28 I hold throughout the exchange, 

Ai(A) = Ao(A) 



(31) 



and the states with rii = and 1 acquire identical Berry 
phases under the exchange. 



This does not imply that the exchange is trivial, since non- 
trivial effects can still arise due to explicit differences be- 
tween the initial and final states. Recall that by assumption 
the Hamiltonian returns back to its original form after the ex- 
change. Therefore we must have /i(Afinai) — e'^^ /i(Ainit) 
and da{Xfinai) — e'^°(iQ,(Ai„it), for some unspecified phase 
factors 9f and 9^- By choosing the bulk quasiparticle wave- 
functions appropriately, we can always set without loss of 
generality 0^ — 0, though we must allow for nontrivial 6f 
since we restricted /i(A) to satisfy Eqs. (28 1 for all A. The 



initial and final ground states are then related by 



|0n2 • • -nAf)/ 
|ln2 • ■■nM)f 



|0n2 • • • nM)i 
l"-2 ■ ■■nM)i- 



(32) 



The most important implication of these results is that the 
outcome of the exchange follows solely from the evolution of 
the localized Majorana wavefunctions for 71.2, which deter- 
mines Of. [This is true up to an overall Berry phase arising 
from Aj(A) which, by contrast, depends additionally on the 
evolution of the bulk quasiparticle states. But this contribution 
is non-universal in any event as we will see later] It is worth 
emphasizing that the wavefunctions corresponding to 71.2 de- 
pend neither on the presence of the additional Majoranas in 
the system, nor on changes to the bulk quasiparticle wave- 
functions far from 71 2. This leads to a great simplification: 
if we know how 71 2 transform in some minimal configura- 
tion, then the same transformation must hold when arbitrarily 
many additional distant Majoranas are introduced by, say, nu- 
cleating other topological regions into the network. Once this 
is known we can determine the unitary operator U12 that acts 
on 71 2 to implement their exchange. Importantly, U12 is ba- 
sis independent as it only acts on the Majorana operators, and 
therefore allows one to deduce the evolution of the many-body 
ground states under the exchange (up to an overall phase) in 
any basis one chooses. 

In the more general situation where one exchanges 7^ and 
7j, exactly the same analysis holds when one works in a ba- 
sis where one of the occupation numbers corresponds to the 
ordinary fermion operator /' = (7^ + i^j)/2, as briefly men- 
tioned at the beginning of this subsection. The implications 
discussed in the preceding paragraph that follow from these 
results also carry over to this more general case. We stress 
that although it is easiest to draw these conclusions by choos- 
ing a different basis depending on which Majoranas are being 
braided, one again ultimately deduces basis-independent op- 
erators Uij that act on 7^ to implement their braid. The set 
of unitary operators {Uij} are sufficient to deduce how the 
degenerate ground states transform in a fixed basis under ar- 
bitrary exchanges (up to an overall non-universal phase); see 
for example Eq. (9) in Ref.il3i 

We have now distilled the problem of demonstrating non- 
Abelian statistics down to that of determining how the Majo- 
ranas 7i and 72 transform under the exchanges shown in Fig. 
|3] These two elementary operations alone indeed allow one to 
braid arbitrary pairs of Majorana fermions in networks com- 
posed of trijunctions, such as that shown in Fig. 4(a) of the 
main text. We proceed now by analyzing these two cases in 
turn. 
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N+1 



2N+1 



N 



FIG. 6; Lattice structure giving rise to the T-junction. 



3. Exchange o/ 71 _ 2 from Figs, )-(d) 

Consider the T-junction in Fig. [3ja). As in the main text 
we will describe each wire by Kitaev's toy lattice model. (We 
remind the reader that the experimentally realistic semicon- 
ductor wire Hamiltonian for the junction can be smoothly de- 
formed to this minimal model without closing the gap. We 
choose to work with the latter since it is more convenient, 
and are free to do so because exchange statistics is a univer- 
sal property that is insensitive to details of the Hamiltonian 
as long as the gap remains finite.) We take the horizontal 
and vertical wires to consist of 2N + 1 and N sites as Fig. 
|6]illustrates, and denote the spinless fermion operators for the 
horizontal and vertical chains respectively by and Cy. For 
simplicity we will take the pairing and nearest-neighbor hop- 
ping amplitudes to be equal, which allows us to conveniently 
express the Hamiltonian as 



H 



2N+1 N 

— ^ y ] MxC^jCj; — y ^ jlyCyCy 

2N 

x=l 
N 



(33) 



where we identify cn+i = cm+i- 

It is important to recall from the discussion below Eq. ( [TT] i 
that the superconducting phases for the T-junction are only 
well-defined with respect to a direction in each wire. We have 
chosen to write the pairing terms above as —t^e'^'^CxCx+i + 
h.c. and similarly for the terms. The former of course equals 
—txe^^'^^^'^ Cx+iCx+h.c, so one can always equivalently refer 
to the superconducting phase in the horizontal wire as (/) with 
respect to the 'right' direction or + tt with respect to the 
'left' direction. Similarly, one can label the superconducting 
phase for the vertical wire as either (f) with respect to the 'up' 
direction or + tt with respect to the 'down' direction. One 
can even divide up each wire into distinct regions and use a 
different convention in each one. We will indeed find this use- 
ful to do here. Throughout this subsection we will employ a 



convention where the superconducting phases are labeled as 
(/) + vr in the left half of the horizontal wire, (f) in the right half, 
and ^ + TT in the vertical wire. This will prove extremely con- 
venient for describing the Majorana operators as we will see 
shortly. Note also that here and in the following subsection 
we will assume that </> and do not differ by an integer mul- 
tiple of TT so that we avoid generating tt junctions during the 
exchange. 

In the initial configuration the vertical wire is non- 
topological, while the horizontal wire is topological and thus 
exhibits end-Majoranas 71 2 which we would like to exchange 
counterclockwise as in Fig. [3|a)-(d). Given our discussion 
in the previous subsection, we will choose 71 2 to satisfy the 
usual relations 7^ = 7j and 7^ = 1 throughout, for in this 
case we need only compare the initial and final Majorana op- 
erators to deduce the outcome of the exchange (up to an over- 
all non-universal phase). In the main text we already deduced 
how the Majoranas transform under this exchange, but here 
we will rederive this transformation rule in the case where the 
superconducting phases are held fixed. In other words, we 
now wish to braid 71,2 by only varying jix, fj-y,tx, and ty. To 
implement the exchange we take these couplings to depend on 
a parameter A that varies from Aa, corresponding to the initial 
setup of Fig. 3[a), to A^, corresponding to the final configura- 
tion of Fig.13 d). The intermediate values At c will correspond 
to Figs.[3jb) and (c). 

We will model the initial setup at A = Aa by assuming for 
simplicity that tx = t > 0, Hx = iy = 0, and jly = ji < Q. 
We can then read off the form of the initial Majorana operators 



fromEq. (33 1: 



71 (AJ = e-*(^+^)/M + /i.e. 



72 (A.) 



-i0/2„t 



1 + h.c. 



(34) 
(35) 



since these combinations are explicitly absent from the initial 
Hamiltonian. Note that the overall signs of 712 (Aa) have been 
chosen completely arbitrarily above. Suppose now that we 
adiabatically transport 71 rightward and downward, leading to 
the configuration of Fig. |3jb). For example, 71 can be trans- 
ported rightward one site by taking /ii — A/iandii — {l — \)t 
and varying A from to 1 . Upon similarly transporting 71 all 
the way to site 1 of the vertical chain, at A = A^ we end up 
with iix — tx = in the non-topological region of the hor- 
izontal wire, iJ^x — and tx — t > in the topological region 
of the horizontal wire, and fly — and ty = t. Once again, 
we can then read off the Majorana operators. 



7i(Ab) 
72(^6) 



72(Aa), 



(36) 
(37) 



where si = ±1. At this point the utility of our convention for 
the superconducting phases becomes clear: with this labeling 



scheme, the phase factors appearing in Eqs. ( 34 1 to ( 37 1 follow 



directly from the superconducting phase 'felt' locally by the 
Majoranas. 

In contrast to the overall signs of 71,2 (Aa) above, si and the 
sign of 72 (Aft) are not arbitrary since these operators evolve 
adiabatically from 7i,2(Aa). Since we did not transport 72 
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in this step, clearly this operator must remain invariant as we 
have written. Determining the sign si is, however, trickier 
One can in principle determine si by tracking the evolution of 
the Majorana wavefunction by explicit calculation. This brute 
force approach does not provide much insight, however, so we 
will deduce the sign si by different means. Specifically, we 
will ask what happens when we go from Fig. [3ja) to (b) in 
a slightly deformed Hamiltonian which smoothly connects to 
the one we discussed above. 

Suppose that instead of the superconducting phases jump- 
ing discontinuously between the horizontal and vertical wires, 
the phase in the latter varied spatially from + tt at the bot- 
tom to + TT at the top. (All the phase variation can happen 
very locally near the trijunction, so long as it is smooth.) Let 
us denote the Majoranas in this deformed problem by ji 2 
which can be expressed as 

7^(A) = e"^'^^^/^Ci{\)^ + h.c. (38) 
72(A) = e-"^^(^)/2^2(A)^ + /i.c., (39) 

where Ci.2(A) correspond to the fermion operators for the 
sites where 2 localized at the 'time' A of the exchange. 
The phases (j)i,2{X) are simply given by the superconducting 
phases felt locally by the Majoranas, and thus rotate as we 
transport 2- Upon going from Fig. l3la) to (b), 02(A) of 
course remains unchanged while 0i(Ajrotates from + tt 
to + TT. Whether this rotation happens clockwise or coun- 
terclockwise, however, makes all the difference: if in the 
former case e-^f-^+^^^Z^ ^ g-i(0+7r)/2^ ^^^^ jj^g j^jjgj. 

g-i(0+7r)/2 _^ 



'i{4>+Tr)/2 



This follows directly from the 
fact that e~*'*j '^'^'/^ changes sign upon shifting (t>j by 27r. Be- 
cause of this mathematical property, it is useful to introduce 
branch cuts in the space of 0j(A), precisely as in Ivanov's 
construction'^*, so that the Majorana operators are single- 
valued away from the cut but change sign when <j)j (A) crosses 
the cut. Arbitrarily, we will take the branch cut to occur at 
0, = 0. 

If one is to avoid generating spurious zero modes at a tt 
junction during this step (and thus smoothly connect to the 
problem of interest), there is a unique orientation in which the 
superconducting phase in the vertical wire can vary. Namely, 
it must vary from + tt at the bottom to </) + tt at the top with- 
out passing through cj). [To understand why going through (p 
would generate a tt junction with zero modes, recall the dis- 



cussion below Eq. (Ill, keeping in mind our labeling conven- 
tions.] Consequently, (/)i(A) rotates from + tt at A = Aq 
to + TT at Afc, also without passing through (p. In Fig. [T^a) 
for example, the rotation happens clockwise as shown by the 
solid red line. Since the deformed problem smoothly connects 
to the physical problem of interest throughout this step of the 
exchange, we can now easily deduce the sign si that deter- 
mines 71 (Af,): si — —1 if 01 (A) crosses the branch cut while 
Si — +1 otherwise. 

Next, let us return to the original problem and similarly 
transport 72 leftward, producing the configuration of Fig.[3jc). 
At A = Ac, we again have fi^ = fJ- < 0, tx — in the 
non-topological region of the horizontal wire, fi^ = and 
tx = t > in the topological region of the horizontal wire. 




(b) (f)2{X) 




FIG. 7: Trajectory of (a) 0i(A) and (b) 02(A) defined in Eqs. 1 38 1 
and {i9\ , which are the superconducting phases 'felt' locally by the 
braided Majoranas 7 j and 72 of the deformed Hamiltonian described 
in the text. The solid red lines correspond to the trajectories obtained 
for the type of exchange illustrated in Figs.[3ja)-(d), while the dashed 
red lines correspond to Figs. |3|e)-(h). Because the Majorana opera- 
tors depend on e"^J we introduce a branch cut along 0j = so 
that 7j 2 are single-valued away from the cut but change sign when 
0j crosses the branch cut. In the exchange from Figs.[3|a)-(d), both 

01 and 02 always rotate by vr with the same orientation, and hence 
only one of the two Majorana operators acquires a minus sign due to 
the branch cut. In the exchange of Figs.[3fe)-(h), however, 0i and 

02 rotate by tt with opposite orientations, so either both 7^ 2 acquire 
a minus sign due to the branch cut or neither do. There is an addi- 
tional minus sign though that 72 picks upon crossing the trijunction 
between Figs.|3[f) and (g), so the net result is that in both exchanges 
one of the two Majoranas changes sign, just like for vortices in a 2D 
p + ip superconductor. Since the deformed model and the Hamilto- 
nian of interest can be smoothly connected without a gap closure, the 
same is true for the original Majoranas, 71,2- 



and fly = and iy = t. The Majorana operators then evolve 



to 



71 (Ac) 
72 (Ac) 



71 (Ab) 



S2 



(40) 
(41) 



for some sign S2- 

To determine S2, we again turn to a deformed problem 
which provides useful insight and eschews the need for brute 
force calculation. As before, we seek a companion Hamil- 



tonian with Majoranas 2 defined as in Eq. (38 1 and (39i, 
where the superconducting phase 02(A) felt locally by 72 
varies smoothly as 72 shuttles leftward. This can be achieved 
by deforming the original Hamiltonian so that the supercon- 
ducting phase in the horizontal wire varies spatially from 0+7r 
on the left, to in the center where it meets the vertical wire. 
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and then to cj) on the right. There is again a unique orienta- 
tion in which this variation can take place if this problem is to 
smoothly connect with the original one — the superconducting 
phase must never pass through (f) + it lo avoid generating a tt 
junction in the topological region during this step. Given how 
the horizontal wire's superconducting phase varies in the de- 
formed Hamiltonian, it follows that (t>2{^) must rotate from 
at Af, to + TT at Ac in an orientation that passes through cj). In 
the example of Fig.|7|b), this implies that 02 rotates clockwise 
as shown by the solid red line, crossing the branch cut (wavy 
line) in the process. We can now trivially extract the sign S2 
using the fact that the deformed Hamiltonian and the physi- 
cal Hamiltonian of interest connect smoothly to one another 
throughout this step: S2 = —1 if 02(A) crosses the branch cut 
but S2 = +1 otherwise. 

Going back to the original problem, let us complete the ex- 
change by transporting 71 up and to the right, producing the 
configuration of Fig.[3|d). The Hamiltonian returns to its orig- 
inal form specified above, and the Majoranas evolve to 



72 (Ad) 



SlSl 

72 (Ad) 



J + h.c. 



(42) 
(43) 



for some sign s'^ . 

One can deduce s'l using the same approach employed to 
find si. We analogously introduce a deformed Hamiltonian 
where the superconducting phase varies spatially in the ver- 
tical wire from + tt at the bottom to (p at the top, without 
passing through + vr to avoid generating zero modes at a tt 
junction. The Majorana j-^ then feels a local superconducting 
phase 01 (A) which rotates from + tt at Ac to at Ad with an 
orientation that avoids passing through + tt. In the example 
of Fig. |7ja), this rotation happens clockwise as indicated by 
the solid red line. Using continuity of the deformed prob- 
lem and the physical problem we immediately deduce that 
s'l = —1 if 01 (A) crosses the branch cut during this step and 
s'l = +1 otherwise. 

At the end of the exchange we find that 71 (Ad) = 
sisi72(Aa) and 72(Ad) = S27i(Aa)- In Fig.|7]only the phase 
02 crosses the branch cut as one can see from the solid red 
lines; in this case si = s'l = 1 and S2 = —1 so that one 
obtains 71 — > —72 and 72 — > 71. More generally, it follows 
from our discussion above that 0i and 02 both always rotate 
by TT with the same orientation, and thus necessarily only one 
of the two crosses the branch cut. Thus we obtain the result 
that the exchange sends 



71 
72 



S72 
-S71 



(44) 
(45) 



with some sign s that depends on 0, 0, and where one chooses 
the branch cut. (We can always get rid of s by absorbing this 
sign into the definition of 72.) Thus we have derived the trans- 
formation rule for the exchange of Fig. |3ja)-(d) in a comple- 
mentary manner to that of the main text. 

It is worth briefly contrasting these two approaches. In this 
subsection we studied the exchange of two Majorana fermions 
when the superconducting phases in the wires were kept fixed. 



We deduced how the Majoranas transformed here by consider- 
ing a companion problem where the superconducting phases 
were allowed to change, but only locally in the vicinity of 
the trijunction. This approach required more extensive for- 
malism but had the virtue of adhering closer to the physical 
situation for experiment. It is also interesting that the pic- 
ture that emerged is remarkably similar to Ivanov's picture for 
exchange of vortices in a 2D p + ip superconductor '^, de- 
spite the absence of vortices in the wire network. In the main 
text we instead deformed the Hamiltonian such that the su- 
perconducting phases in each segment of the network varied 
globally during the exchange so as to keep the Hamiltonian 
real. At the end, however, the superconducting order param- 
eter changed sign, so to return the Hamiltonian back to its 
original form we reversed the sign of the pairing by sending 
— (71 ~ *72)/2 — > ip ■ This method had the virtue that all 
Berry phases explicitly vanished while the Majoranas were 
transported, allowing one to get to the answer very directly 
and with minimal formalism. We stress that these two ap- 
proaches are not unrelated, and it is interesting to note the con- 
nection between the deformed companion problem introduced 
in this section and the deformed Hamiltonian considered in the 
main text. In the former case the superconducting phases felt 
by the Majoranas ji 2 varied smoothly as they moved across 
the network, and as a result the phases 0i 2 smoothly rotated 
by TT during the course of the braid. The latter case effectively 
corresponds to an extreme version of this wherein 0i 2 acquire 
the entire tt rotation only at the very end of the exchange. 



4. Exchange of'y\^2from Figs.^e)-(h) 

We will now analyze the counterclockwise exchange of 
Figs. [3]^e)-(h), where 71 and 72 reside on different topologi- 
cal regions of the network. The same general approach ap- 
plied in the previous subsection will be applied here. As be- 
fore, the T-junction will be described by H in Eq. ( 33 1, the 



Majoranas 712 will be defined so that 7^ = 7^ and 7I = 1 
throughout, and their exchange will be implemented by tak- 
ing the couplings /i^;, ij.y,tx, and ty dependent on a parame- 
ter A. We will hold and fixed, but will again deduce the 
evolution of the Majoranas by considering a deformed com- 
panion Hamiltonian where the superconducting phases vary 
locally near the trijunction. It will be convenient to employ a 
slightly different labeling convention for the superconducting 
phases here compared to the previous subsection: we will la- 
bel these by in the left half of the horizontal wire, + tt in 
the right half, and + tt in the vertical wire. The 'time' A will 
now vary from Ac to A/i, respectively corresponding to the se- 
tups of Figs. [3];e) and (h); the intermediate values A^ and Xg 
will similarly correspond to Figs.|3jf) and (g). For simplicity, 
as above we evolve the Hamiltonian such that the couplings 
at each of these 'times' are given hy ^.j.,,jiy — ji < and 
ixjiy — in the non-topological regions, while iJ.x,fJ-y — 
and ta:,ty = t > in the topological regions of the network. 

Suppose that at Ag the Majoranas 71 and 72 reside at sites 
X and X' of the horizontal wire, respectively; the operators 
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can then be written 

7l(Ae) 
72 (Ae) 



6-^^/24 + h.C. 

g-»(0+.)/2^t^^ + h.C. 



(46) 
(47) 



[We can safely ignore the other two Majoranas in Fig. |3je), 
since they do not evolve at all during the exchange of 71,2 and 
thus 'factor out' ; see again Sec. IV B 2 for a rigorous discus- 
sion. We can even include a coupling between the left and 
right ends of the horizontal wire if desired to fuse these addi- 
tional Majoranas.] When A increases to A/, we arrive at the 
setup of Fig.[3|f) and the Majoranas evolve to 



71 (A/) 

72 (A/) 



72 (Ae), 



(48) 
(49) 



for some sign a^. Notice that with our modified labeling 
scheme, the phase factors appearing in the above operators 
again simply follow from the superconducting phases felt lo- 
cally by each Majorana. 

One can deduce cri by the usual procedure applied in the 
previous subsection. We introduce a deformed Hamiltonian 
where the superconducting phase in the vertical wire varies 
spatially from (p+n at the bottom to cf) at the top, without pass- 
ing through + TT to avoid generating spurious zero modes at 
a TT junction during this step. We define the Majoranas 2 
for this deformed Hamiltonian as in Eqs. ( [38] l and ( [39] i; the 
phases 01.2 (A) are again precisely the smoothly varying su- 
perconducting phases felt locally by 2- I" particular, as 
7]^ shuttles rightward and then downward, it feels a supercon- 
ducting phase (A) that rotates from i/i at Ae to ^ + tt at A/ in 
an orientation that avoids passing through (/> + tt. In the exam- 
ple from Fig. [vja), (pi rotates from cj) to (p + n counterclock- 
wise, as the dashed red line indicates. Because the deformed 
Hamiltonian and our original Hamiltonian can be smoothly 
connected throughout this step without closing a gap, we can 
immediately deduce the sign (Ti: (Ti = — 1 if 0i crosses the 
branch cut while tJi = +1 otherwise. 

Returning to the original Hamiltonian, suppose we now 
transport 72 leftward, generating the configuration of Fig|3jg). 
The Majoranas then take the form 



7i(A3) = 7i(A/) 



72 (Ag) =0-2 6 



.-i0/2„t 



(50) 
(51) 



for some sign 0-2 • Deducing (T2 is trickier than the other signs 
because midway through this step 72 sits exactly at the trijunc- 
tion, with all three emanating wire segments being topologi- 
cal. To analyze this step we begin by introducing a deformed 
Hamiltonian where the superconducting phase in the horizon- 
tal wire varies spatially from from on the left, to + tt in the 
center, to + tt on the right. To avoid generating unwanted 
zero modes during this step, we require that the superconduct- 
ing phase varies in this fashion without passing through 4>. As 
the Majorana 72 moves leftward and approaches the junction. 



this operator evolves according to Eq. ( 39 1 with the local su- 
perconducting phase 02 (A) varying from + tt to + tt with- 
out crossing 0. 



Eventually 72 sits exactly at the trijunction, and here we 
need to proceed with care. Let us denote this point by A = 
Atri- First, we note that in our deformed Hamiltonian the hor- 
izontal wire forms a tt junction at Atri- Because the vertical 
wire is also topological, however, no additional zero modes 
are generated as discussed around Eq. ([12]). As we also dis- 
cussed there, because of the tt junction the wavefunction for 
72 at Atri has no weight on the vertical wire. We explicitly 
find that here 73 evolves to 



72 (Atri) 



1 



-j(0+ir)/2 



t ) 
N+2) 



h.C. 



(52) 



(This can be deduced by considering only the sites N,N + 
1, + 2 of the horizontal wire and site A^ of the vertical wire.) 
The minus sign appearing in front of cjy is absolutely crucial. 
Because of this extra minus sign, as 72 moves off of the tri- 
junction and proceeds leftward, it will subsequently evolve 
according to 



72(A > Atri) = -e-^*^(^)/2^2(A)^ + h.C. 



(53) 



where 02(A) corresponds to the superconducting phase felt 
locally by 73 during the latter half of this step. Thus 73 feels 
a local superconducting phase 02(A) which varies smoothly 
from 0+7r at A/ to at Xg in an orientation that passes through 
+ TT, but additionally picks up an extra minus upon crossing 
the trijunction. In the example from Fig.jTJ^b), 02(A) rotates 
clockwise as shown by the dashed red line. Using continuity 
between the deformed and original Hamiltonians, we can now 
conclude that (T2 = +1 if 02 crosses the branch cut during this 
step while a2 ~ —1 otherwise. 

Finally, let us return to the original Hamiltonian and com- 
plete the exchange by transporting 71 up and to the right. Thus 
we arrive at the setup of Fig. |3jh), and the Majoranas become 



7i(Ah) 
72(Ah) 



CiCTi 

72(Ag) 



g-»(0+.)/2^t^^ + h.C. 



(54) 
(55) 



for some sign a'l that we can determine by the usual means. 
Introduce a deformed Hamiltonian where the superconducting 
phase in the vertical wire varies from + tt at the bottom to 
+ TT at the top without passing through to avoid additional 
zero modes. As the Majorana for this companion problem 
shuttles up and to the right, it feels a local superconducting 
phase 01 (A) which varies from + tt at A^ to + tt at A/j 
with precisely this orientation. In the example of Fig.jTJa), 0i 
rotates from 0+7r to 0+7r counterclockwise, as the dashed red 
line shows. By continuity, we immediately find that a'l = —1 
if 01 crosses the branch cut during this step while a[ ~ +1 
otherwise. 

The final and initial Majorana operators are related by 
7i(Ah) = cr'iCrij2{K) and72(A,,) = cr27i(Ae). For the ex- 
ample shown in Fig. |7]we obtain (Ti = a[ = +1 and (72 = ^ 1, 
so here 71 — >■ —72 and 72 — > 71. More generally, both 0i and 
02 always rotate by tt with opposite orientations for this type 
of exchange, and so either 71 and 72 both pick up a minus 
sign because of the branch cuts or neither do. The Majorana 
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72 always acquires an additional minus sign, however, upon 
crossing the branch cut, so we find as before that 



7i 

72 



Sl2 

-s-fi- 



(56) 
(57) 



The sign s again depends on 0, and how one orients the 
branch cuts, but is the same for the exchange of Figs. [3ja)- 
(d) and that of Figs. [3|e)-(h). Thus as one would intuitively 
expect, both kinds of counterclockwise braids transform the 
Majoranas in an identical fashion, despite the fact that in one 
case the Majoranas are initially bridged by a topological re- 
gion while in the other they initially reside on disconnected 
topological segments. This proves our claim in the main text 
that these exchanges are equivalent. 



5. Implications fornon-Abelian statistics 

Consider now a network composed of some arbitrary ar- 
rangement of trijunctions, such as that of Fig. 4(a) from the 
main text. The elementary braids of Fig. |3]constitute the basic 
operations needed to exchange Majorana fermions in this set- 
ting. Putting together all the results obtained so far in this sec- 
tion, exchanging counterclockwise any two Majorana modes 
7i and 7^ without disturbing any other Majoranas in the sys- 
tem (apart from, perhaps, trivial translations) sends 



7j -Sijji 

Ik Ik {k^ 



(58) 
(59) 
(60) 



The sign depends on the initial sign choices for the oper- 
ators 7i j, the superconducting phases in the system, and how 
one orients the branch cuts as discussed above. Once these 
are specified, can be deduced for any exchange using the 
simple rules outlined in the previous two subsections. If Uij 
denotes the operator that implements this counterclockwise 
braid, we have 



a. 



V2 ■ 



(61) 



[One can easily apply the analysis of the previous subsections 
to the clockwise analogue of the exchanges of Figs. [3|a)-(d) 
and Figs. |3je)-(h). This of course leads to the result that the 
clockwise braid of 7^ and jj is generated by U^j ~ ^Ij-^ 
Non-Abelian statistics for the wire network now follows from 
the fact that 



(62) 



for i ^ k. 

While trijunctions alone are sufficient to allow for non- 
Abelian statistics, we note in passing that it is of course pos- 
sible to consider more general networks featuring some ar- 
bitrary number of wires meeting at a junction. The general 



results established in Sec. IV B 2 still apply here, though ad- 
ditional cases can arise beyond those considered in Fig. [3] As 



an example, suppose one fabricated a 'H-' junction where four 
wire segments meet at a point. If the entire junction is topo- 
logical, then four Majoranas will generically appear. If we 
exchange a given pair, which of the braided Majoranas ac- 
quires a minus sign can not be immediately deduced from our 
results above, though this case can be analyzed exactly along 
the lines of how we studied the exchanges of Fig. [3] Our aim 
is not to be completely exhaustive here, however, so we do not 
pursue such cases further in this work. 



6. Many-body Berry phase calculation for a system with two 
Majorana fermions I 

Although we have already established non-Abelian statis- 
tics in wire networks, in the final parts of this section we will 
explicitly analyze the evolution of the full many-body ground 
states under exchange in some tractable cases. This will serve 
to not only support our previous analysis, but also enable us 
to discuss important issues such as the overall Berry phase 
acquired by the ground states upon braiding Majoranas. We 
begin here with the simplest case and return to the initial setup 
in Fig. [3]; a). To exchange the Majoranas 712 as in Fig.jSj^a)- 
(d), here we will follow the strategy adopted in the main text 
and keep the Hamiltonian purely real during this exchange 
(until the very end, when we will allow the Hamiltonian to 
become complex). Again, this assumption has the virtue that 
the wavefanctions can then also be chosen real, so that in spite 
of their complex evolution the Berry phase accumulated as the 
Majoranas are transported vanishes identically. 



As in Sec. IV B 3 we will describe the T-junction by the lat- 



tice Hamiltonian in Eq. (33 1. We will revert back for the rest 



of the Supplementary Material to the labeling scheme where 
the superconducting phase is denoted by in the horizontal 
wire and in the vertical wire (with respect to the 'right' and 
'up' directions, respectively). For convenience we will de- 
form the Hamiltonian describing the initial setup of Fig. [3|a) 
to the following: 



N 



2N 



= 1 



"x+1 



),(63) 



with fj, < and t > 0. Here we have set the horizontal wire's 
superconducting phase to (j) — and chemical potential to 
/ij = 0, and turned off the pairing and hopping in the vertical 
wire. The Hamiltonian then exhibits only real matrix elements 
as desired. We graphically denote the initial superconducting 
phase in the horizontal wire by the rightward-pointing arrow 
in Fig.[3|a) (a leftward-pointing arrow would indicate a phase 
of TT, which would also keep the Hamiltonian purely real). 

The first term in Hi implies that all Cy fermions in the ver- 
tical wire will be absent in the initial ground states, while the 
second can be recognized as Kitaev's toy model in the special 
limit where n — 0, t — \A\. The end Majorana fermions for 
the horizontal wire take on a particularly simple form in this 
limit, allowing the initial wavefunctions to be easily obtained. 
To do this, we follow Kitaev^^ and decompose in terms of 
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Majorana fermions ^a/b,x via 
1 



B,. 



(64) 



which allows the Hamiltonian to be written as 

N 2N 

-P^'^^l^y - ^^'^lB,xlA,x+i- (65) 



x = l 



The zero-energy end Majorana fermions ja,i and jb,2N+i 
which do not appear in H can be combined into an ordinary 
zero-energy fermion 



while the gapped bulk states are captured by operators 
1 



:{lA,x+i + ilB.x)- 



In terms of d-r, Hj becomes 



N 2N 

H, = -fiJ2clcy+tJ2i'2didx-l) 

y=l x=l 



(66) 



(67) 



(68) 



The end Majoranas give rise to two degenerate initial 
ground states whose evolution we are interested in: |0) i which 
dond annihilates and — dl^^\0)i. The former can be writ- 
ten |0)i — dcnd n^=i dx\ya.c), where |vac) denotes the vac- 
uum of Cj. and Cy fermions. After some algebra, the normal- 
ized ground states can be written explicitly as 



|0), 



^ N 2N+1 

E E 



p=0 il<---<i2p + l 



■c-Jvac) 



1 

2^ 



N 2N+1 



1+E E 

p=l ii <---<i2p 



|vac). (69) 



Note that we have multiplied |0)i and by overall phase 
factors to make each wavefunction purely real. Although the 
ground states have different fermion parity, both yield the 
same average particle number 



N = 



2iV + 1 



(70) 



corresponding to half-filling of the horizontal chain. 

Let us now transport the Majorana fermions as outlined 
in Figs. |3ja)-(d), keeping the Hamiltonian (and ground state 
wavefunctions) real and avoiding spurious zero-energy along 
the way. For example, 71 can be transported rightward one 
site by adding the following term to Hi, 



8H 



A/icjci - Xt{c\ + ci)(4 - C2) 



(71) 



(with < 0) and varying A from to 1. As usual, as we 
so transport 71 and 72 we must avoid having two neighboring 



topological regions whose superconducting phases differ by 
TT, for in this case a pair of 'accidental' zero-energy Majorana 
modes appears at the junction. It is therefore useful to employ 
arrows as shown in Figs. |3ja)-(d) to signify the sign of the 
pairing in each topological region. Two inward or two out- 
ward arrows meeting at the junction correspond to a tt junc- 
tion and must be avoided. Figures [3|a)-(d) illustrate that in 
accordance with this simple rule, we can indeed swap the po- 
sitions of 71 and 72 while keeping the Hamiltonian and wave- 
functions purely real, consequently acquiring no Berry phases 
whatsoever in the process. However, the arrows and hence the 
sign of the pairing in the topological region unavoidably re- 
verse, as seen by comparing Figs.|3ja) and (d). Thus we have 
not yet completed an exchange in the usual sense. 

At this stage we have adiabatically evolved the Hamiltonian 

to 



JV 



2N 



H' = -/i^£t£^-t^(4-c,)(4+i+c,+i),(72) 



x=l 



corresponding to Hi with the sign of the pairing reversed, and 
the wavefunctions to 



JV 2^+1 

|o)' = 2i^E E ("ircU---< 

p=0 ii<---<i2p+l 

N 2N+1 



vac) 



1 

2^ 



1+E E (-iK 

p=l ii<---<i2p 



•4 

I2p «1 



|vac(J73) 



[Modulo phase factors, these wavefunctions can be obtained 
by sending Cx icx in Eqs. ([69|.] To complete the exchange, 
let us now return the Hamiltonian to its original form by adia- 
batically rotating the superconducting phase in the topological 
region from tt back to 0. The Hamiltonian then involves com- 
plex matrix elements, which implies that Berry phases need 
no longer vanish here. As we will see, however, the Berry 
phase contributions for this final step can be easily calculated. 
To this end, consider 



N 2N 
y—1 x—1 



J\it/2 



Cx) 



(74) 



Upon varying A from to 1, the superconducting phase rotates 
by TT such that H{\ = 0) = H' and H{\ = 1) = Hi as 
desired. The ground states of H{X) are 



10(A)) = 



^E E (-ire-'-(-+^ 



/2) 



p=0 ii<---<i2p4 



X C, 



i2p+i 

xe 



11(A)) 



X c. 



2N 



■ |vac; 



N 2N+1 



1 + E (-l)^e~*^^f 

p=l ii<---<i2p 



I vac). 



(75) 
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Importantly, |0(A = 0)) = |0)' and |1(A = 0)) = so that 
the wavefunctions evolve smoothly throughout. Note also that 
we have inserted an arbitrary phase factor 9 above. We will 
select this phase momentarily such that the Berry phase ac- 
quired by each wavefunction during this final stage also van- 
ishes. The outcome of the exchange is then simpler to inter- 
pret, since one simply compares the initial states |0)i and \ 
with the final states |0>/ = |0(A 1)> and = |1(A = 1)). 

Using Eqs. ( [75] l, one can now compute the Berry phases; 
we find 

Im / dA(0(A)|aA|0(A)) = Im / dA(l(A)|aA|l(A)) 
Jo Jo 

2 



(76) 



(Off-diagonal components such as (0(A) |9a| 1(A)) vanish triv- 
ially due to the different fermion parity of the ground states.) 
This result is quite sensible given that both wavefunctions de- 
scribe on average N/2 Cooper pairs whose phase rotates by 
IT. We now choose 



Ntt 



(77) 



so that the Berry phases vanish as desired. Only the explicit 
relative phases between the initial and final wavefunctions re- 
main. For A = 1 the factors of (—1)^ cancel in Eqs. (75 1, 
yielding 



|0>/ = 



|0). 



(78) 



Crucially, the ground state |1) acquires an additional phase 
factor of i relative to |0) under the exchange. Neglecting an 
overall phase factor, the unitary operator that generates this 
relative phase can be written 



Ul2 



72 71 



(79) 



where we have identified 71 = ja.i and 72 = Jb,2N+i- This 
coincides with the expression obtained in the main text by 
somewhat different means, and is identical to the unitary op- 
erator generating the exchange of vortices in a spinless p + ip 
superconductor' 

Three important comments are warranted here. First, it is 
worth emphasizing again that in practice one need not perform 
any rotation of the superconducting phases to exchange Majo- 
ranas or realize non-Abelian statistics, as we have already seen 
in the preceding subsections. We analyzed the problem in this 
way solely because the many-body wavefunctions and Berry 
phases could be computed very easUy in this approach. In the 
physical situation appropriate for quantum wires, the effective 
p-wave superconducting phases in the horizontal and vertical 
wires will differ by tt /2, and to implement the exchange one 
only needs to apply local gate voltages to exchange the Ma- 
joranas. One can also compute the Berry phases in this situa- 
tion (which we will indeed do momentarily in a simple case), 
though the calculations are much more complicated. 



Second, while a specific overall phase has been calculated 
in Eqs. ( [78] l, this phase is certainly non-universal. It depends 
both on the precise way in which one exchanges the Majo- 
ranas as well details of the Hamiltonian, and is therefore not 
particularly meaningful in this context. For instance, had we 
rotated the superconducting phase from tt back to before 
moving 71 all the way to the right in Fig.[3|d), a different over- 
all phase would emerge. (As we found earlier, the average 
number of particles encoded in the wavefunctions when the 
rotation takes place affects the Berry phase and would be dif- 
ferent in this case.) Furthermore, as we demonstrate below the 
overall phases depend on the specific form of the Hamiltonian 
even when the superconducting phases remain fixed; see Eqs. 
( [ST] ) and ( |105| l. This result is perhaps not too surprising — the 
overall phases follow from many -body wavefunctions that en- 
code not only topological information, but also non-universal, 
high-energy physics. 



Third, to obtain the result in Eq. (79 1 we chose in Eq. (74i 
to rotate the superconducting phase counterclockwise from tt 
to in the final step of the exchange. Had we alternatively 
chosen to rotate the phase in a clockwise fashion, the ground 
state |1) would pick up a relative phase of —i instead of +i 
under the exchange compared to |0). Interesting physics un- 
derlies the ambiguity. To understand this, first note that the 
cases where one rotates the superconducting phase from tt to 
counterclockwise versus clockwise differ by an overall ro- 
tation by 2tt. In any superconductor, rotation of the super- 
conducting phase by 2tt effectively changes the sign of all 
fermion operators; in particular, here such a 2tt rotation sends 
71.2 ^71, 2- Remarkably, exchanging 71 and 72 twice 
(with the same orientation) while keeping the superconduct- 
ing phases fixed also sends 712 — > —71.2- In other words, 
braiding 71 all the way around 72 is equivalent — modulo an 
overall phase — to a braidless operation wherein the positions 
of the Majoranas remain fixed but the superconducting phase 
advances by 2tt. Consequently, if under a counterclockwise 
exchange |1) picks up a relative phase of +i compared to |0), 
then subsequently rotating the superconducting phase by 2tt 
effectively converts this into a clockwise exchange with |1) 
acquiring a relative phase of —i compared to |0) . 



We should emphasize that precisely the same conclusions 
apply to a 2D spinless p + ip superconductor featuring two 
vortices. In that case, however, this observation is less inter- 
esting. With more than two vortices, rotating the overall phase 
of the 2D p + ip order parameter by 2tt changes the sign of 
all the Majoranas, which is not very useful. In wire networks, 
however, by fabricating a series of Josephson junctions along 
the wires one can in principle wind the superconducting phase 
only along the region supporting a particular pair of Majo- 
ranas, thus effectively implementing pairwise braids in a po- 
tentially useful — though not topologically protected-way. We 
also note that such "braidless exchanges" were recently dis- 
cussed in a rather different setting by Teo and Kane^^. 
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7. Many-body Berry phase calculation for a system with two 
Majorana fermions II 

If one transports the Majoranas while leaving the supercon- 
ducting phases fixed (as one would do in practice), no such 
ambiguity arises. To demonstrate this we now consider a four- 
site model, with the geometry of Fig.[6]with iV = 1, described 
by the following Hamiltonian: 

H = -nic\ci - fisclcs - fiic\ci 

+ ti{c\ + Ci)(4 - C2) +^2(4 + C2)(4 - C3) 



+ ?i(e' 



T/4J 



C2)(e 



5i).(80) 



As before, ci.2,3 correspond to sites on the horizontal chain, 
while ci couples to C2, forming the vertical bond of the T- 
junction. The ti_2 and ti terms represent nearest-neighbor 
tunneling and pairing with equal amplitude, while iii 3,fli 
denote the chemical potentials (we will never need a chem- 
ical potential for C2 and so have excluded such a term above). 
Importantly, the superconducting phases have been set to — 
and will remain fixed at — zero on the horizontal bonds and 
—tt/2 on the vertical bond, similar to the physical situation 
for semiconducting wires as discussed in the main text. 
We wish to evolve the coupling constants in Eq. (8O1 to 



implement an exchange of two Majorana modes as shown in 
Figs. [3|a)-(d). It is simplest to carry out the exchange piece- 
wise in six stages. We construct the many-body ground states 
such that they are continuous between steps and compute the 
Berry phases separately within each step. Although in prin- 
ciple the wavefunctions can be obtained analytically, the ex- 
pressions are algebraically very complicated and not terribly 
illuminating. Thus we will only outline the calculation here. 

In the first step, the couplings evolve according to ti = 
A's = 0, ^2 = = — Ml = — A/i, and ti = t(l — A), 

with t, n > and A varying from to 1. In the initial configu- 
ration with A = 0, Majorana fermions 71 and 72 respectively 
reside on sites 1 and 3 of the horizontal chain, realizing a con- 
figuration analogous to Fig.jsj^a). We define / = (71 +1^2)/'^ 
such that the initial ground state |0)i annihilated by / has odd 
fermion parity, while the even-parity initial ground state is 
= /^|0)i. Ramping up A — >■ 1 shuttles 71 rightward to 
site 2. Since = throughout this step, the Hamiltonian is 
real and the wavefunctions can thus also be chosen real. The 
Berry phase therefore vanishes trivially here. 

In the second step, the couplings evolve according to ti = 
^3 = 0, ^2 = t, /ii — jli — —(1 — A)/i, and ti — Xt. 
Varying A from to 1 now shuttles 71 downward, realizing a 
configuration analogous to Fig. [3|b). In the third and fourth 
steps we similarly evolve the couplings to transport 72 from 
site 3 to site 2, and then from site 2 to site 1, leading to the 
configuration of Fig.[3|c). The fifth step transports 71 upward 
to site 2. Note that in steps 2 through 5, the Hamiltonian ex- 
hibits complex matrix elements, and non-zero Berry phases 
thus emerge. Finally, in step 6 71 moves rightward to site 
3, completing the exchange. In the last step the Hamiltonian 
can once again be chosen real, but we explicitly introduce A- 
dependent phase factors into the ground states [as we did in 



eels the Berry phase contributions from steps 2 through 5. All 
the physics is then contained in the relative phases between 
the initial and final wavefunctions; we obtain 



|0>/ = -^e-''-/m^ 



\1) 



f 



-i37r/4| 



1) 



(81) 



Unlike the case above where we rotated the superconducting 
phases during the exchange, here the relative phase factors 
have been determined unambiguously. Most importantly, the 
ground state |1) picks up an extra phase factor of i compared 
to |0) as expected. 



8. Evolution of many-body ground states under exchange in a toy 
model exhibiting four Majorana modes 

Finally, we will examine how the many-body ground states 
transform under an exchange of the type outlined in Figs.[3]^e)- 
(h), where there are four Majoranas present. We would in par- 
ticular like to address here the following puzzle. Suppose we 
proceed as we did earlier for the exchange of Figs. [3|a)-(d) 
and transport 71 2 while keeping the Hamiltonian purely real 
to avoid Berry phase accumulation. It is once again instructive 
to view the sign of the pairing in the topological regions with 
arrows as displayed in Figs.|8ja)-(d). As the figure illustrates, 
it is now possible to exchange 71 and 72 without reversing the 
sign of the pairing in the process. In other words, we can keep 
the Hamiltonian purely real, swap the locations of the Majo- 
ranas, and return the Hamiltonian back to its original form — 
without closing a gap. One might therefore worry that this 
type of exchange is trivial, but this is not so. This can be eas- 
ily deduced from the perspective of the Majorana operators 
using the preceding results from this section: upon crossing 
the trij unction 72 acquires an additional minus sign, and so 
under the exchange we get the usual result that 71 — > 572 
and 72 — S71, for some sign s. From the perspective of the 
wavefunctions, how then do the ground states possibly acquire 
the relative phase factor of i that follows from this transforma- 
tion? 

The subtlety arises because there are now four Majorana 
fermions rather than two, and in this case one gets less mileage 
out of keeping the Hamiltonian real during the exchange. To 
illustrate the point, let 71.2 denote the Majoranas we wish to 
exchange, and 73 4 the stationary Majoranas of Figs.|3je)-(h). 
Defining 



fA 



(71 + ^72) 

(73 +i74), 



(82) 



Eqs. (75 1] such that the Berry phase from step 6 exactly can- 



we see that there are now two degenerate ground states in each 
fermion parity sector: |0, 0) which both Ja and /s annihilate, 
=/i/l|0,0), |1,0) =/l|0,0),and|0,l) =/l|0,0). 
Reality of the Hamiltonian does not imply that these four 
ground states can each be chosen real. Indeed, we provide 
an example below where this is clearly not possible. Rather, 
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71 
(a) 



72 



72 





FIG. 8: Exchange of two Majorana fermions separated by a non- 
topological region. Here we envision transporting tlie Majoranas 
while keeping the Hamiltonian purely real. The superconducting 
phase in the topological regions can then only take on two values, 
or TT, which we indicate by arrows above. Unlike in the exchange of 
Figs.|3[a)-(d), here we can exchange Majorana fermions while keep- 
ing the Hamiltonian purely real, maintaining the gap, and returning 
the Hamiltonian back to its original form (i.e., without reversing the 
sign of the pairing). As explained in the text, this does not mean that 
the exchange is trivial; indeed, the Majoranas transform just as they 
do in the braid of Figs.|3[a)-(d). 



this condition only guarantees reality of specific linear com- 
binations of these ground states, which in general can vary as 
the exchange takes place. In other words, the reality condition 
does not preclude the phases of the above ground states from 
evolving nontrivially during the exchange (see below for an 
explicit example). Drawing conclusions about the exchange 
from this route therefore requires a more detailed analysis than 
in the case with only two Majorana fermions. To remedy this 
issue, one might be tempted to modify the setup of Figs.[3je)- 
(h) by connecting the horizontal wire into a loop, then fusing 
73 and 74 so that only the two Majoranas which we exchange 
remain. One will quickly discover, however, that in this case 
the positions of 71 and 72 can not be swapped while keeping 
the Hamiltonian real and all other excitations gapped. Specif- 
ically, in the process one necessarily ends up with a config- 
uration similar to Fig. |8jc), except with two arrows pointing 
either into or out of the junction; that is, one can not avoid tt 
junctions here. 

Let us consider now a toy problem that provides an illustra- 
tive minimal setting in which the wavefunctions during such 
an exchange can be analyzed explicitly. Specifically, we ex- 
amine the four-site setup shown in Fig.|9]and described by the 
following purely real Hamiltonian; 

H = -Ai4c4+i(cl +ci)[(4-C2) + (4-C3)] 

+ [i24(4+C2)+<34(4+C3)](4-C4), (83) 

with t > 0, fi < 0, and t24,i34 > 0. In spite of the small 
number of sites, this Hamiltonian supports four zero-energy 
Majorana modes for any values of /i, ^24, and t34. To get in- 
tuition here, it is useful to think of site 1 as forming a tt junc- 
tion between sites 2 and 3. This gives rise to two Majorana 



(a) 1 



72 72 



^1 • 



72 72 



(d) 1 

'29 '2. ^9 



71 



FIG. 9: Minimal four-site setup that supports four Majorana modes 
(only 7i and 72 are shown for clarity). The Hamiltonian is chosen 
so as to exchange 71 and 72 as sketched in (a)-(d), mimicking the 
exchange of Fig. [3|e)-(h) in a tractable setup. Solid lines denote 
bonds with non-zero pairing whose sign is indicated by the arrows. 



modes which are independent of the parameters appearing in 

(84) 



Eq. (83 1. One of these, 73, resides at site 1; 

73 = «(4 - Ci). 
The other, 74, resides on sites 2 and 3; 



74 



V2 



(4 



C2 - cl + C3 



(85) 



The locations of the second pair of Majoranas, 71 and 72, de- 
pend on /i, i24, and i34. We will vary these parameters so 
as to carry out an exchange of 71 and 72 in a manner that is 
analogous to the exchange of Figs.[3|e)-(h). 

To help establish a connection between the setup of Fig. 
[3je) and the present four site-problem, imagine first forming 
a loop out of the horizontal wire so that the two topological 
regions connect. The outer Majoranas of Fig. [3]^e) — which 
are analogous to 73.4 in our setup — can then be generated by 
forming a tt junction in the topological region. The four-site 
problem shrinks this tt junction to the smallest possible size. 
The other two Majoranas 71 2 will initially be separated by an 
unpaired region, similar to the non-topological segments con- 
necting 712 in Fig.[3]^e). We carry out their exchange piece- 
wise in three stages to reduce the algebraic complexity of the 
problem. Care will be taken to ensure that the wavefunctions 
and operators defined below evolve continuously in between 
each of these stages. 

(1) In the first stage, we evolve the Hamiltonian by taking 



^ = (1 - A)^4 

^24 — At 
^34 = 



(86) 



(with /i4 < 0) and varying A from to 1 . Initially when A = 0, 
71 and 72 are situated at sites 2 and 3, respectively, as Fig.|9|a) 
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illustrates. Ramping up A to 1 shuttles 71 to site 4, leading to 
the configuration of Fig.[9];b). More precisely, 71 2 are given 
by 

-fl = — [-A*(4+C2)+2t24(4+C4)] (87) 
024 



-3-rC3, (88) 

where we have defined a24 = ^^^^24: + M^- The finite-energy 
fermion operators which annihilate the ground states are 

4 = ^[^/2(cI+ci) + (4-C2) + (4-C3)] (89) 

= [2*24(4 + C2) + + "24)4 + (a^ - "24)04]- 

2a24 

(90) 

Now define and analogously to Eqs. (82 1. Suppress- 



ing the A dependence for notational simplicity, the four de- 
generate ground states are then 



|0,0)/ 



(91) 



|i,o)/ = /rio,o)/ 

|0,1)/ = fl\0,0)i, 

where |0, 0)/ is annihilated by fs, d^, and dg. With the 
above definitions and some time to carry out the algebra, one 
can obtain these ground states for arbitrary A. When A = 
leading to the initial configuration shown in Fig. |9|a), the 
wavefunctions are 
1, 



|0,0). 

|1,0). 
|0, 1). 



'44 +44] I vac) 



+ 44 



^44 



34] 

-e'^c\ - ic\ - 4 - e-*^44cl] 



(92) 



vac) 



"4 + j444]|vac) 



Clearly none of these can be made real by introducing overall 
phase factors (though one can readily verify that a purely real 
basis does exist by considering linear combinations of these 
states). When A ^ 1 and we arrive at the configuration shown 
in Fig.|9jb), the wavefunctions evolve to 

|0,0)fc = ^ 



2V2 

+ 44' 



'44 + s 



f44 



e *'icjjC344]l'v^'^) 



|l,l)h = 



2^/2 



+ 44 



^44 



'^44 + e*^44 + i4444]|vac) (93) 



(II) For the second stage of the exchange, we evolve the 
Hamiltonian according to 



^24 
*34 





{l-\)t 

xt. 



(94) 



Here varying A from to 1 leaves 71 unchanged but adiabati- 
cally transportS72 from site 3 to site 2, leading to the config- 
uration of Fig.bfc). Defining /3 = \/t^4^ + t\^, the Majorana 
fermion operators at this stage obey 



C4 



-^34(4 + C2) + ^24(4 + C3)], 



while the gapped quasiparticle operators are 



d-B = ^[i24(4 + C2)+<34(4 + C3)+/3(4 



(95) 
(96) 

(97) 

C4)r98) 



(Note that there is some freedom for how one implements this 
step. The expression for 7!^ does not depend on the specific 
parametrization of t24 and ^34 above. If desired, one can first 
turn on ^34 while keeping ^24 fixed, and then turn off <24- The 
end result in either case is the same.) The wavefunctions can 
again be obtained for arbitrary A after some tedious algebra. 
In particular, when A — > 1 bringing the system to the setup of 
Fig.|9|c), the wavefunctions evolve to 



1 



2^-1 



e 44 + s 



+ c'^cl, + 44 ^ *44 ^ e** 4444]|vac) 



44 

11,1). = -^\e 



''4'--3^2'--lJ 



2^2 



+ i44 + 44 + 44 + 



+ e"*T44 + e*^44 + j4444]|vac) (99) 



|l,0)c = 



1 



|0,1). 



2^2 
1 



^4 + 4 



6*2444 + 6 



5444 +444] |vac) 



*444 

Notice how the phase factors in the wavefunctions evolve non- 



'^444]|vac) 



trivially in passing from Eqs. (93 1 to (99 1, despite the reality 
of the Hamiltonian. 

(Ill) To conclude the exchange, we now choose 



|l,0)fc = 



2\/2 

6*^444 



-6**4 - ic\ 



e '4444 



-3^2"-l 



'^444 



■444] I vac) 



|0,l)h = 



2V2 



c +e 



-if 



ci +e 



-if 



+ *444 



- 444 



*444 



6*^0444] I vac) 



M 

*24 
*34 



A^4 



(l-A)t 



(100) 



and again vary A from to 1. In this final step, 72 remains 
unchanged while 71 moves adiabatically to site 3 as in Fig. 
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|9|d). The Majorana operators now obey 



— [-M(4+C3) + 2*34(4 + C4)] (101) 



"34 



72 



(102) 



with — \/4:t^4 + /i^, and the gapped quasiparticle opera- 
tors are 



jiii 



(103) 



2a 



— [2*34(4 + C3) + (^ + a34)4 + (Ai - "34)24] 
134 



(104) 



Computing the ground states as before, we obtain the expected 
result that the final and initial ground states are related by 



|0,0)/ 

|1,1>/ 
|1,0)/ 
|0,1>/ 



|0,0)^ 
z|l,0), 

|o,i).. 



(105) 



That is, the ground states with an Ja = (71 + *72)/2 fermion 
present acquire a phase factor of i under the exchange. Up 
to an overall non-universal phase factor, this transformation is 
generated by the unitary operator U12 = e'^''^'^'-^^, as obtained 
earlier 

To close this section, we remark that one may object that in 
obtaining this result we have simply compared the initial and 
final states. Since the above wavefunctions are not real, one 
may in particular ask whether the exchange is tainted by Berry 
phases. It is not — it is always possible to simply change to a 
real basis by suitably superposing these wavefunctions, and 
in such a basis the absence of Berry phases is manifest. The 
exchange indeed is governed solely by the difference between 
initial and final states. 



Derivation of the fractional Joseplison effect in a simple 
model 



For pedagogical purposes, we will review here the 'frac- 
tional Josephson effect' originally predicted by Kitaev^^ and 
discussed by other authors in the context of ID wires' ^ and 
other topological systems'^ We will examine this ef- 

fect in a minimal setup where all calculations can be explicitly 
carried out, although the qualitative aspects of the physics are 
more universal. Consider two topological superconducting 
wires forming a Josephson junction as shown schematically 
in Fig. 



10 a). The phases of the p-wave order parameters are 



taken to be 



in the left/right wires, which are coupled by 



a weak (compared to the gap in each wire) electron tunneling 
term at the junction. The full Hamiltonian reads 



H = Hr.+H„ + Hr 



(106) 



where /r describe the left/right regions and Hr represents 
the electron tunneling term coupling the wires. For compu- 
tational simplicity, we model the left and right regions as A^- 
site chains described by Kitaev's toy model with fi = and 



(a) 




(b) 
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FIG. 10: (a) Schematic of the Josephson junction formed by two 
topological wires with p-wave superconducting phases 4>l/r- The 
wires couple at the junction through an electron tunneling term with 
strength T. (b) Bogoliubov-de Gennes spectrum as a function of 
{(t)L — ipR) /{2n) for the effective Hamiltonian in Eq. l |l 17| cho- 
sen to describe the junction. Only the solid lines denote physically 
distinct states. The states centered around E/{2t) = 1 represent 
ordinary bulk quasiparticles, while the state near zero energy repre- 
sents the quasiparticle formed when the two end Majoranas at the 
junction fuse. The energy and hence Josephson current correspond- 
ing to the latter exhibit 4n periodicity in — (I>r- The ordinary 
bulk quasiparticle states, however, contribute only to the usual 2-k- 
periodic Josephson effect. 



t = |A|. In this case we have 



x=l 



X (e" 



-#o/2„t 



-^a.x+1 



i0„/2 N 
^ ^a,2;-|-l J ; 



(107) 



(108) 



for a — L/R, along with a tunneling term 

Hr = -r(4,ArC_R,i + h.c). 

When F = 0, two Majorana fermions reside at the junction; 
turning on F 7^ generally fuses these to an ordinary finite- 
energy quasiparticle state. We wish to compute the zero-bias 
current flowing across the junction, 

I ^-Y^^'^iN^R^i+h-c), (109) 

in the ground state as well as the excited state where this quasi- 
particle state is occupied. 

We proceed by first diagonalizing Hj^ in the usual way. 
Writing c^^.x = e-"^°/2(-^a^ _,_ ^7^^)72 and then defining 
da,x = hXx+i + *7s,2;)/2. one obtains 



Hn 



N-1 



(24 



1) 



(110) 



22 



It is useful to group the end Majorana fermions residing at the 
junction into an ordinary fermion operator via 



1 



4nd = o(7a,i +*7s.jv)- 



The tunneling term, which we will treat as a perturbation, can 
then be written 



r 



Hr = -{C[d^j^^^{dl^^_,+dL,N-i) + h.c.] 



7t 



+ <S'Mlnd('^kl ~ '^^'l + d,L,N-l + dL,N-l) + h.c] 
+ C(2dLd4nd - 1)} (112) 



with 



C ^ cos(A0/2) 
5 = sin(A(/)/2) 

A(j) = (l)L- (jjR- 



(113) 
(114) 
(115) 



Rewriting the expression for the current in this basis, one ob- 
tains the familiar relation 



2e d{Hr) 
h dAcj) 



(116) 



Notice that the fermion operators (iL,i...7v-2 and dB..2...N-i 
essentially drop out from the problem — the full Hamiltonian 
separately conserves the fermion number for each of these 
states and they do not contribute to the Josephson current. 
Thus for the purposes of evaluating the current, the prob- 
lem maps onto a simpler Hamiltonian involving only d^nd, 
dL,N-i, and du,!- In terms of dA — {dL,N-i + di?,i)/v^ 
and ds = {dL^^-i — dR^i)/^/2, this effective Hamiltonian 
becomes 



i/cff = t[{2d\dA - 1) + {2dldB - 1)] + Hr, (117) 
where now 

Hv = ^{C[{2dl^^d,^d-l) + {d\dA-dldB)] 

+ C(44 + h.c.) + V2S[idl^^{d\ + de) + h.c.]}. 

(118) 

Applying degenerate perturbation theory to obtain the ener- 
gies of the dend, dA, and fermions to 0[(r/t)-^], we obtain 



7t 



HcS 



EA(f\fA-l]+EB (flfB-l: 



^end [ /cjid-^cnd (-> 



(119) 



The operators /a/b /end correspond to states that evolve from 
d-A/B /end due to the tunneling perturbation. Their energies to 
the desired order are 

Ea/b = 2t±^cos(A0/2) + ^[5-3cosA(/)j;i2O) 



We can now evaluate the Josephson current in the ground 
state, as well as the excited state where the finite-energy quasi- 
particle formed from the fused Majoranas is occupied. Equa- 



(111) tions ( 116 1 and ( 1 19 1, along with the above energies, yield 



= ±|sin(A0/2)-^sinA0, 



(122) 



where the +/— sign corresponds to the current obtained when 
the /ond fermion is unoccupied/occupied. The second term 
represents the standard Josephson current that is 2tt periodic 
in Acj). This contribution reflects Cooper-pair tunneling and 
thus arises at second-order in perturbation theory. More inter- 
estingly, the first term exhibits Att periodicity and has a topo- 
logical origin since it arises solely from the Majoranas fused 
at the junction. This contribution reflects ?i first-order process 
corresponding to single-electron tunneling, which is possible 
at zero bias because the Majoranas form a zero-energy state at 
the junction when F = 0. 



It is interesting to observe from Eq. (120i that the Ja/b 
fermions also pick up a first order correction to their energy 
from r. Thus one can view each of these states as individ- 
ually contributing both 27r- and 47r -periodic Josephson cur- 
rents. Their 47r-periodic contributions exactly cancel one an- 
other, however, so that only the fused Majoranas contribute 
to this effect. Mathematically, this can be understood from 
the particle/hole-symmetric spectrum of the Bogoliubov-de 
Gennes Hamiltonian in Eq. (|1 This is plotted versus 



A0/ (27r) in Fig. 10 'b). Here only the solid lines denote phys- 
ically distinct states since those with energy E and —E are not 
independent. As one turns on the tunneling strength F from 0, 
the ordinary fermionic states that begin at energy 2t split with 
opposite sign at first order, and only yield a net change in en- 
ergy at second order in F/t. Thus they contribute only to the 
usual 27r-periodic Josephson current. The state of affairs for 
the end state which begins at zero energy is very different — its 
energy also shifts at first order, but its 'partner' which shifts 
in the opposite direction does not represent a physically dis- 
tinct state. It therefore produces an observable 47r-periodic 
Josephson current. 

Finally, it is useful to ask whether the crossings in the spec- 
trum of Fig. 10 b) at A(j) = tt are stable. In the case of the or- 



£;end = Tcos{A(t>/2). 



(121) 



dinary /a,b quasiparticle states, they are certainly not. For in- 
stance, adding a weak superconducting pairing between jv 
and Cfl 1 at the junction lifts the crossings near E — ±2t in the 
figure. When this happens, there is no sense in which the bulk 
quasiparticle states even individually contribute a 47r-periodic 
current. The crossing at £' = 0, however, is stable provided 
the Majoranas at the outer ends of the wires do not overlap 
with those at the junction^^. (The location of the crossing^ 
though need not occur exactly at A^ — tt). This can be un- 
derstood as follows. As long the as occupation number of the 
fermion corresponding to the outer Majorana end states re- 
mains fixed, the ground states at Acj) = and A(p = 2tt have 
different fermion parity. In the former case f^^-^^f end = in 
the ground state while in the latter /^jj^/end = 1- If this cross- 
ing could be removed, then one would be able to adiabatically 
evolve A(f> from to 27r while remaining in the ground state, 
but this can not happen unless the outer end Majoranas trans- 
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fer a fermion to the junction. It is useful to keep this in mind 
when considering the left path in Fig. 4(b) of the main text, 
where 72 crosses the junction. At fixed phase A(^* ^ tt, the 
ground state is always accessible here, precisely because 72 
overlaps with the fused Majoranas 734 at the junction during 
this process. However, when 72 resides at the far right end of 
the wire, then (neglecting residual overlap between the Ma- 
joranas) the ground state will no longer accessible when the 
phase difference changes by 27r. This explains why the frac- 
tional Josephson current in Eq. (9) of the main text exhibits 
27r periodicity in the initial phase difference Aip^, but Air pe- 
riodicity in A (p. 
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